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ABSTRAGT 

The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty 
star-forming galaxies (DSEGs) over a wide range in redshift. A detailed understanding of these DSEGs 
is hampered by the limited spatial resolution of Herschel. We present 870 /im 0'/45 resolution imaging 
from the Atacama Large Millimeter/submillimeter Array (ALMA) of 29 HerMES DSEGs with far- 
infrared (EIR) flux densities in between the brightest of sources found by Herschel and fainter DSEGs 
found in ground-based sub-millimeter (sub-mm) surveys. We identify 62 sources down to the 5cr 
point-source sensitivity limit in our ALMA sample (a ^ 0.2 mJy), of which 6 are strongly lensed 
(showing multiple images) and 36 experience significant amplification (/i > 1.1). To characterize the 
properties of the ALMA sources, we introduce and make use of uvmcmcfit, a publicly available 
Markov chain Monte Garlo analysis tool for interferometric observations of lensed galaxies. Our lens 
models tentatively favor intrinsic number counts for DSEGs with a steep fall off above 8mJy at 
880 jam. Nearly 70% of the Herschel sources comprise multiple ALMA counterparts, consistent with 
previous research indicating that the multiplicity rate is high in bright sub-mm sources. Our ALMA 
sources are located significantly closer to each other than expected based on results from theoretical 
models as well as fainter DSEGs identified in the LADOGA EGDES Submillimeter Survey. The high 
multiplicity rate and low projected separations argue in favor of interactions and mergers driving the 
prodigious emission from the brightest DSEGs as well as the sharp downturn above Sggo = 8 mJy. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: high-redshift 
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1. INTRODUCTION 

Galaxies selected in blind surveys at far-infrared (FIR) 
or sub-millimeter (sub-mm) wavelengths are generally 
known as dusty star-forming galaxies (DSFGs; for a re¬ 
cent review, see Gasey et al. 2014). They cover a wide 
range in redshift from z^0.5toz>6 (Chapman et al. 
2005; Casey et al. 2012a; Messias et al. 2014; Riechers 
et al. 2013), with a significant component at z ^ 2 (Casey 
et al. 2012b; Bothwell et al. 2013), when they represent 
the most FIR-luminous objects in existence during this 
epoch. They are usually signposts of significant over¬ 
densities (Daddi et al. 2009; Capak et al. 2011) (c.f. Rob¬ 
son et al. 2014) and likely represent the formative stages 
of the most massive elliptical galaxies found in the lo¬ 
cal Universe (e.g., Ivison et al. 2013; Fu et al. 2013). 
Moreover, they constitute an important component of 
the overall galaxy population at z ^ 2 (e.g., Magnelli 
et al. 2011), when the star-formation rate density in the 
Universe peaked (e.g., Lilly et al. 1996; Madau et al. 
1996). 

Our collective understanding of DSFGs is currently 
taking a dramatic leap forward thanks in large part to 
the Herschel Space Observatory {Herschel; Pilbratt et al. 
2010). Herschel has revolutized the size and depth of 
blind surveys at FIR wavelengths. In particular, the iifer- 
5c/ie/Multi-tiered Extragalactic Survey (HerMES; Oliver 
et al. 2012) and the Herschel Astrophysical Terahertz 
Large Area Survey (H-ATLAS; Eales et al. 2010) to¬ 
gether have surveyed ^ 650 deg^ at 250 /rm, 350 /rm, and 
500 jam to the confusion limit of Herschel (a ~ 6 — 7 mJy 
in each band Nguyen et al. 2010), plus an additional 
^ 350 deg^ to a shallower level (approximately double the 
confusion limit). A similar effort to survey large areas of 
the sky has been undertaken at longer wavelengths by the 
South Pole Telescope (SPT; Carlstrom et al. 2011) and 
the Atacama Cosmology Telescope (Swetz et al. 2011). 

Theoretical expectations based on the redshift distri¬ 
bution and luminosity function of DSEGs suggested that 
HerMES and H-ATLAS would be efficient tools for dis¬ 
covering strongly lensed DSEGs (e.g., Blain 1996; Ne- 
grello et al. 2007). Submillimeter Array (SMA; Ho et al. 
2004) imaging at 870 jam with sub-arcsecond resolution 
has confirmed this, with > 85% of the brightest sources 
found by Herschel that satisfy S'soo > 100 mJy being 
gravitationally lensed by an intervening galaxy or group 
of galaxies along the line of sight (Negrello et al. 2010; 
Conley et al. 2011; Riechers et al. 2011a; Bussmann 
et al. 2012; Wardlow et al. 2013; Bussmann et al. 2013). 
Sources discovered in SPT surveys have also been shown 
to have a high probability of being strongly lensed (Vieira 
et al. 2013; Hezaveh et al. 2013). However, statistical 
models significantly over-predict the median magnifica¬ 
tion factor experienced by a Herschel DSEG of a given 
S'soo (Bussmann et al. 2013). This could herald new 
insights in our understanding of the bright end of the 
intrinsic DSEG number counts or in the nature of the 
deflectors. 

We here present Atacama Large Millime- 
ter/submillimeter Array (ALMA) Cycle 0 imaging 
at 870/im of a sample of 29 HerMES DSEGs. Three 
aspects of our dataset make it uniquely suited to improv¬ 
ing our understanding of the bright end of the intrinsic 
DSEG number counts. Eirst, the sample occupies a 


distinct regime in flux density between the brightest 
Herschel DSEGs (almost all of which are lensed) and 
much fainter DSFGs found in ground-based surveys 
(most of which are expected to be unlensed; e.g., Hodge 
et al. 2013). Second, the ALMA images are extremely 
sensitive (rms point source sensitivity of a ~ 0.2 mJy) 
and all 29 HerMES DSFGs are detected (which was 
not the case in previous similar studies with shallower 
imaging; e.g., Smolcic et al. 2012; Barger et al. 2012; 
Hodge et al. 2013). Third, the typical angular resolution 
is 0745 and nearly all sources detected by ALMA are 
spatially resolved. 

We also obtained Gemini-South optical imaging 
to complement our existing set of ancillary multi¬ 
wavelength imaging. We use those data in this paper to 
identify lensing galaxies, which are typically early-types 
with little on-going star-formation and therefore exhibit 
very weak sub-mm emission. 

In Section 2, we characterize our sample and present 
our ALMA and Gemini-South imaging. Section 3 
presents our model fitting methodology and model fits 
for all ALMA sources (lensed and unlensed) using UVM- 
CMCFIT, a publicly available modified version of the 
visibility plane lens modeling software used in Bussmann 
et al. (2012, 2013). Results on the effect of lensing for the 
observed properties of the Herschel DSFGs in our sam¬ 
ple, as well as the multiplicity rate and typical angular 
separation between sources after delensing the ALMA 
sources, appear in Section 4. We scrutinize statistical 
predictions for the magnification factor at 870 jam {jisio) 
as a function of the fiux density at 870 jam (Ssjo) and dis¬ 
cuss implications for the bright end of the DSEG number 
counts in Section 5. Finally, we present our conclusions 
in Section 6. 

Throughout this paper, we assume a flat cosmology 
with Hq =69 km s“^ Mpc“^, = 0.29 (Hinshaw et al. 
2013). 

2. DATA 

In this section, we describe the selection of our Herschel 
DSEG sample, present our ALMA high-spatial resolution 
imaging of thermal dust emission, and present Gemini-S 
optical imaging that we use to identify intervening galax¬ 
ies along the line of sight. 

2.1. Selection of DSFG Sample 

The starting point for the sample selection is source 
extraction and photometry. For the objects in this pa¬ 
per, individual catalogs were generated for each of the 
250 /im, 350 /im, and 500 jam Herschel Spectral and Pho¬ 
tometric Imaging REceiver (SPIRE; Griffin et al. 2010) 
channels using the SUSSEXtractor peak finder algorithm 
(Savage & Oliver 2007). Our sample includes 29 DSFGs 
drawn from five independent, confusion-limited fields 
in HerMES with declinations below +2° and totaling 
55deg^. 

The sample was selected to be the 29 brightest DS¬ 
FGs in the Southern sky that are not known radio AGN, 
nearby late-type galaxies, or Galactic emission. The se¬ 
lection was designed to assemble a large sample of lensed 
galaxies in the ALMA-accessible HerMES fields, and was 

https://github.com/sbussmann/uvmcmcfit 
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constructed from the SUSSEXtractor catalogs (Smith 
et al. 2012), which were available prior to the ALMA Cy¬ 
cle 0 deadline. Subsequently, improved efforts to deblend 
SPIRE photometry at 500 /im using StarEinder (Wang 
et al. 2014) were introduced that formed the basis of the 
lens selection criteria used in Wardlow et al. (2013). As 
a result of the improved deblending algorithms in the 
StarEinder catalogues and Wardlow et al. (2013) a num¬ 
ber of objects in our sample have significantly lower S '500 
values in the StarEinder catalog than in the original SUS- 
SEXtractor catalogs. This and further investigation into 
the StarEinder catalogs shows that their original flux was 
boosted by blending with nearby sources rather than by 
gravitational leasing. Eor this reason, the objects in this 
sample comprise a combination of lenses and blends of 
multiple sources. 

We used positional priors based on the ALMA data 
presented in this paper to obtain the best possible es¬ 
timates of the total SPIRE flux densities for each Her¬ 
schel source. We also used Spitzer/MIPS (Rieke et al. 
2004) imaging to take into consideration the presence of 
nearby 24/rm sources that are not detected by ALMA 
but may still contribute to the 250 jam emission detected 
by Herschel. Additional details on our methodology are 
provided in Appendix A. The SPIRE flux densities mea¬ 
sured in this way represent our “fiducial” flux densities 
and are presented in Table 1. Interested readers may re¬ 
fer to Table 5 for a comparison of the fiducial, StarEinder 
and SUSSEXtractor flux densities in tabular form. 

Eigure 1 shows that the Herschel-ALMA sample is set 
clearly apart from the very bright Herschel DSEGs that 
are selected to have S '500 >100 mJy and have been shown 
to be almost entirely lensed DSEGs (Negrello et al. 2010; 
Wardlow et al. 2013; Bussmann et al. 2013). In contrast, 
the sample in this paper is expected to include a mix 
of lensed and unlensed DSEGs. On the other hand, the 
HerMES survey area is 200 times larger than that of the 
Large Apex Bolometer Camera Extended Chandra Deep 
Eield Survey (LESS WeiB et al. 2009). This explains why 
the median S'soo in our sample is ^ 4 times brighter than 
the median S '500 in the sample of ALMA-detected sources 
in LESS, known as ALESS (Swinbank et al. 2014). Our 
Herschel-ALMA sample opens a new window of discovery 
space on the bright end of the DSEG number counts. 

In detail, two of the sources in the Herschel-ALMA 
sample (HXMMOl and HXMM02) overlap with the “con¬ 
firmed lensed” sample in Wardlow et al. (2013) as well 
as with the Herschel-^MA sample in Bussmann et al. 
(2013). A further eight appear in the “Supplementary 
sample” of Wardlow et al. (2013). The remainder have 
5^500 < 80 niJy and thus do not appear in Wardlow et al. 
(2013). 

Table 1 provides reference data for the Herschel-ALMA 
sample, including centroid positions measured from the 
ALMA 870 yum imaging (see Section 2.2). The centroid 
positions serve as the reference point for subsequent off¬ 
set positions of lenses and sources described in later ta¬ 
bles. This is a useful choice (rather than the SPIRE cen¬ 
troid or ALMA phase center) because it minimizes the 
number of pixels needed to generate a simulated model 
of the source and therefore minimized memory and cpu 
usage when lens modeling. 

2.2. ALMA Observations 
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Figure 1. Herschel/SPIKE photometry of all galaxies in the Her¬ 
MES phase I catalog with declination < -h2 ° and signal to noise ra¬ 
tio greater than 5 at 350 jim and 500 /j,m (log of number of galaxies 
shown in grayscale). The sample of HerMES sources in this paper 
are shown with green squares (“Hersc/ie/-ALMA”). The very bright 
Herschel DSFGs from Bussmann et al. (2013) (“Hersc/ie/-SMA”) 
are shown by red circles, and lensed SMGs discovered by the SPT 
that have published lens models (“SPT-ALMA”) are represented 
by cyan stars (Hezaveh et al. 2013). A magenta diamond shows the 
location in this diagram of the stacked signal from ALESS DSFGs. 
Representative error bars are shown in the lower left corner. The 
Herschel-ALMA sample fills the gap in 500/rm flux density space 
between 50-100 mJy. 


ALMA data were obtained during Cycle 0 over 
the from 2012 June to 2012 December (Program 
2011.0.00539.S; PI: D. Riechers). The observations were 
carried out in good 870 jam weather conditions, which 
resulted in typical system temperatures of Tgys ~ 130 K 
and phase fluctuations of ^ 10°. Each target was ob¬ 
served until an rms point-source noise level near the 
phase center of a 0.2 mJy per beam was achieved. This 
typically required 10 minutes of on-source integration 
time. Eor the observations targeting the CDES, ELAISS, 
and COSMOS fields, the data reach a ^ 0.14 mJy per 
beam. The number of antennas used varied from 15 to 
25. The antennas were configured with baseline lengths 
of 20 m to 400 m, providing a synthesized beamsize of 
^ 0'.'5 X 074 EWHM while ensuring that no flux was re¬ 
solved out by the interferometer (since our targets all 
have size scales smaller than 1" — 2". When possible, 
track-sharing of multiple targets in a single track was 
used to optimize the uv coverage. 

The quasars J0403—360, J2258—279, B0851+202, and 
J2258—279 were used for bandpass and pointing calibra¬ 
tion. The quasars J0403—360, J0106—405, J0519—454, 
J1008+063, and J0217-I-017 were used for amplitude and 
phase gain calibration. The following solar system ob¬ 
jects were used for absolute flux calibration: Callisto 
(CDES targets); Neptune (XMM targets); Titan (COS¬ 
MOS targets); and Uranus (ADES and XMM targets). 
Eor HELAISS02, no solar system object was observed. 
Instead, J2258—279 was used for absolute flux calibra¬ 
tion, with the flux fixed according to a measurement 
made two days prior to the observations of HELAISS02. 

All observations were conducted with the correlator 
in “Erequency Domain Mode”, providing a total usable 
bandwidth of 7.5 GHz with spectral windows centered at 
335.995 GHz, 337.995 GHz, 345.995 GHz, 347.996 GHz. 
We searched for evidence of serendipitous spectral lines 
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but found none (typical sensitivity is a ~ 8mJybeam“^ 
in 15 kmsec”^ bins). Given that our observations cover a 
total of 217.5 GHz in bandwidth, the lack of lines seems 
more likely to be due to limited sensitivity than limited 
bandwidth. 

We used the Gommon Astronomy Software Applica¬ 
tions (GASA, version 4.2.1) package to re-reduce the 
data provided by the North American ALMA Science 
Genter (NAASG). We found that the quality of the pro¬ 
cessed data from the NAASG was very high. However, 
we achieved a significant improvement in the case of the 
ADFS and XMM targets by excluding data sets with 
moderate Tgys and poor phase fluctuations. For a handful 
of targets with peak signal-to-noise ratio (S/N) greater 
than 20, we obtained a ~ 10% improvement in S/N by 
using the GASA selfcal task with the clean component 
model as input to improve the phase gain corrections. Fi¬ 
nally, we updated the absolute flux calibration to use the 
Butler-JPL-Horizons 2012 solar system models 

For imaging, we used the GASA Clean task with 
Briggs weighting and “robust = H-0.5” to achieve an op¬ 
timal balance between sensitivity and spatial resolution. 
We selected the multi-frequency synthesis option to op¬ 
timize uv coverage. We designed custom masks for each 
target in GASA to ensure that only regions with high 
S/N were considered during the cleaning process. 

Figure 2 presents our ALMA images (color scale) in 
comparison to the Herschel SPIRE images (black-white 
contours) originally used to select the targets and noted 
in each panel as either 250 /im, 350 /im, or 500 /im. Each 
panel is centered on the phase center of the ALMA ob¬ 
servations of that target and a white circle traces the 
EWHM of the primary beam of an ALMA 12 m antenna 
at 870 yum. All flux density measurements given in this 
paper have been corrected for the primary beam by divid¬ 
ing the total flux density by the primary beam correction 
factor at the center of the source. This is a valid approach 
because all sources have sizes < 1", such that the vari¬ 
ation in the primary beam correction factor across the 
source is insignificant. A white dashed box represents 
the region of each image that is shown in greater detail 
in Eigure 3. 

In most targets, the peak of the SPIRE map is spa¬ 
tially coincident with the location of the ALMA sources. 
In one case where two ALMA sources are separated by 
^ 10" (HADES08), the elongation in the SPIRE 250/im 


map is consistent with the angular separation of the two 
ALMA counterparts. Otherwise, the SPIRE imaging is 
consistent with a single component located at the cen¬ 
troid of the ALMA sources. This result is not a sur¬ 
prise, given the typical angular separation of the ALMA 
sources (< 5") and the EWHM of the SPIRE beam 
at 250/im (18.1"). We identify and catalog by-eye all 
sources with peak flux density greater than 5cr. 


2.3. Gemini-South Imaging 

Optical imaging observations using the Gemini 
Multi-Object Spectrograph-South (GMOS-S; Hook 
et al. 2004) were conducted in queue mode during the 
2013B semester as part of program GS-2013B-Q-77 (PI: 
R. S. Bussmann). The goal of the program is to use 
shallow ii, r, i, and z imaging to identify structure 
at redshifts below unity and determine which of the 
ALMA sources are affected by gravitational leasing. 
Nearly half of the ALMA sources he in regions with 
existing deep optical imaging, thanks to the extensive 
multi-wavelength dataset available in the HerMES fields 
— these were excluded from our Gemini-S program. 
The remaining targets are: HADES03, HADES08, 
HADES09, HADESIO, HADES02, HADES04, 
HADESOl, HADESll, HELAISS02, HXMMll, 
HXMM12, HXMM22, HXMM07, HXMM30, and 
HXMM04. Each of these targets were observed for a 
total of 9 minutes of on-source integration time in each 
of 14, r, 1 , and z. The observations were obtained 
during dark time in adequate seeing conditions (image 
quality in the 85th percentile, corresponding to ~ 1.1"). 

The data were reduced using the standard IRAE Gem¬ 
ini GMOS reduction routines, following the standard 
GMOS-S reduction steps in the example taken from the 
Gemini observatory webpage 

We used the Sloan Digital Sky Survey (SDSS) or the 2 
Micron All Sky Survey (2MASS) to align the Gemini-S 
images to a common astrometric frame of reference. This 
imposes an rms uncertainty in the absolute astrometry 
of 072 and 074 for SDSS and 2MASS, respectively. The 
astrometrically calibrated Gemini-S images served as the 
basis for aligning higher resolution, smaller field-of-view 
imaging from HST or Keck (when available), which were 
originally presented in Galanog et al. (2014). 


Table 1 

Observed positions and flux densities of ALMA sources. For each Herschel 
source, we give the flducial flux density in all SPIRE bands (see main text), 
as well as the observed positions and flux densities of all ALMA sources. The 
statistical rms and synthesized beam EWHM in each ALMA map (ctalma 
and f^ALMAj respectively) are also given. Positional uncertainties (for 
unlensed sources) range from 07005 for well-detected sources to 07l5 for 
the faintest soures in our sample. Uncertainties in flux density do not include 
the absolute calibration uncertainty of 10 %. Quoted uncertainties in 
Herschel photometry are dominated by confusion noise (Nguyen et al. 2010). 


lAU address^ Short name 

RAsto 

(J2000) 

Decgyo 

(J2000) 

*8250 

(mJy) 

8'350 

(mJy) 

*8500 

(mJy) 

*8870 

(mJy) 

f^ALMA 

(mJy) 

^ALMA 
(" X ") 

Lens 

grade^ 

J003823.6-433707 HELAISS02 

00:38:23.587 

-43:37:04.15 

115 ±6 

124 ±6 

108 ±6 

20.11 ±0.45 

0.14 

0.54 X 0.44 

— 

— SourceO 

00:38:23.762 

-43:37:06.10 

— 

— 

— 

9.22 ±0.17 

— 

— 

C 

— Sourcel 

00:38:23.482 

-43:37:05.56 

— 

— 

— 

4.34 ±0.16 

— 

— 

C 

— Source 2 

00:38:23.313 

-43:36:58.97 

— 

— 

— 

4.16 ±0.32 

— 

— 

c 


28https://sciencemrao.edu/facilities/alma/aboutALMA/Technology/Al2Wy4,e;yy^^ie|/^lma|j)4^|a^g^|^^^^^_^^^_^gg^j^g/p^^^ggg;^g_ 

software/getting-started#gmos 
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Table 1 — Continued 


lAU address^ 

Short name 

RAgyo 

(J2000) 

Dec870 

(J2000) 

‘S '250 

(mjy) 

‘S '350 

(mJy) 

‘S '500 

(mJy) 

(mJy) 

C^ALMA 

(mJy) 

^ALMA 
(" X ") 

Lens 

grade^ 

— 

Source3 

00:38:23.803 

-43:37:10.46 

— 

— 

— 

2.40 + 0.19 

— 

— 

C 

J021830.5-053124 

HXMM02 

02:18:30.673 

-05:31:31.75 

78 + 7 

122 + 8 

99 + 7 

63.33 + 0.58 

0.20 

0.49 X 0.37 

A 

J021841.5-035002 

HXMM31 

02:18:41.613 

-03:50:03.70 

102 + 6 

94 + 6 

65 + 6 

10.80 + 0.46 

0.20 

0.49 X 0.37 

— 

— 

SourceO 

02:18:41.520 

-03:50:04.72 

— 

— 

— 

6.79 + 0.37 

— 

— 

C 

— 

Sourcel 

02:18:41.700 

-03:50:02.57 

— 

— 

— 

4.01 + 0.26 

— 

— 

c 

J021853.1-063325 

HXMM29 

02:18:53.111 

-06:33:24.65 

97 + 6 

102 + 6 

78 + 6 

7.28 + 0.45 

0.20 

0.49 X 0.37 

— 

— 

SourceO 

02:18:53.118 

-06:33:24.19 

— 

— 

— 

5.46 + 0.30 

— 

— 

c 

— 

Sourcel 

02:18:53.095 

-06:33:25.21 

— 

— 

— 

1.82 + 0.38 

— 

— 

c 

J021918.4-031051 

HXMM07 

02:19:18.417 

-03:10:51.35 

89 + 7 

107 + 8 

85 + 7 

29.16 + 0.58 

0.21 

0.49 X 0.38 

A 

J021942.7-052436 

HXMM20 

02:19:42.783 

-05:24:34.84 

72 + 6 

85 + 6 

66 + 6 

17.49 + 0.74 

0.20 

0.49 X 0.37 

— 

— 

SourceO 

02:19:42.629 

-05:24:37.11 

— 

— 

— 

7.15 + 0.44 

— 

— 

X 

— 

Sourcel 

02:19:42.838 

-05:24:35.11 

— 

— 

— 

3.52 + 0.41 

— 

— 

X 

— 

Source2 

02:19:42.769 

-05:24:36.48 

— 

— 

— 

3.42 + 0.26 

— 

— 

X 

— 

Source3 

02:19:42.682 

-05:24:36.82 

— 

— 

— 

2.46 + 0.47 

— 

— 

X 

— 

Source4 

02:19:42.955 

-05:24:32.22 

— 

— 

— 

0.94 + 0.18 

— 

— 

X 

J022016.5-060143 

HXMMOl 

02:20:16.609 

-06:01:43.18 

179 + 7 

188 + 8 

134 + 7 

29.56 + 0.46 

0.20 

0.48 X 0.37 

— 

— 

SourceO 

02:20:16.648 

-06:01:41.93 

— 

— 

— 

16.13 + 0.31 

— 

— 

c 

— 

Sourcel 

02:20:16.571 

-06:01:44.56 

— 

— 

— 

11.56 + 0.32 

— 

— 

c 

— 

Source2 

02:20:16.609 

-06:01:40.72 

— 

— 

— 

1.87 + 0.26 

— 

— 

c 

J022021.7-015328 

HXMM04 

02:20:21.756 

-01:53:30.92 

162 + 7 

157 + 8 

125 + 11 

20.03 + 0.47 

0.23 

0.53 X 0.38 

c 

J022029.2-064845 

HXMM09 

02:20:29.140 

-06:48:46.49 

129 + 7 

118 + 8 

85 + 7 

15.30 + 0.36 

0.20 

0.49 X 0.37 

— 

— 

SourceO 

02:20:29.195 

-06:48:48.02 

— 

— 

— 

8.93 + 0.30 

— 

— 

c 

— 

Sourcel 

02:20:29.079 

-06:48:44.86 

— 

— 

— 

6.37 + 0.19 

— 

— 

c 

J022135.1-062617 

HXMM03 

02:21:34.891 

-06:26:17.87 

114 + 7 

134 + 8 

116 + 7 

22.65 + 0.36 

0.21 

0.48 X 0.38 

— 

— 

Sourcel 

02:21:35.124 

-06:26:16.62 

— 

— 

— 

18.42 + 0.36 

— 

— 

c 

— 

Source2 

02:21:35.132 

-06:26:18.02 

— 

— 

— 

2.19 + 0.20 

— 

— 

c 

— 

SourceO 

02:21:35.136 

-06:26:17.28 

— 

— 

— 

2.03 + 0.18 

— 

— 

c 

J022201.6-033340 

HXMMll 

02:22:01.616 

-03:33:41.40 

101 + 7 

104 + 8 

73 + 7 

11.72 + 0.49 

0.20 

0.52 X 0.38 

— 

— 

SourceO 

02:22:01.592 

-03:33:39.42 

— 

— 

— 

8.17 + 0.32 

— 

— 

c 

— 

Sourcel 

02:22:01.629 

-03:33:43.58 

— 

— 

— 

3.54 + 0.36 

— 

— 

c 

J022205.4-070728 

HXMM23 

02:22:05.362 

-07:07:28.10 

128 + 6 

105 + 6 

68 + 6 

2.93 + 0.15 

0.20 

0.48 X 0.37 

X 

J022250.5-032410 

HXMM22 

02:22:50.573 

-03:24:12.35 

101 + 6 

85 + 6 

61 + 6 

10.19 + 0.28 

0.20 

0.49 X 0.38 

c 

J022547.8-041750 

HXMM05 

02:25:47.942 

-04:17:50.80 

103 + 7 

118 + 8 

97 + 7 

17.96 + 0.43 

0.20 

0.50 X 0.37 

c 

J022944.7-034110 

HXMM30 

02:29:44.740 

-03:41:09.57 

86 + 6 

97 + 6 

75 + 6 

22.76 + 0.28 

0.23 

0.50 X 0.38 

A 

J023006.0-034152 

HXMM12 

02:30:05.950 

-03:41:53.07 

98 + 7 

106 + 8 

82 + 7 

15.56 + 0.37 

0.20 

0.50 X 0.38 

c 

J032752.0-290908 

HECDFS12 

03:27:52.011 

-29:09:10.40 

61 + 7 

82 + 6 

81 + 6 

38.78 + 0.56 

0.15 

0.43 X 0.35 

— 

— 

SourceO 

03:27:52.002 

-29:09:12.07 

— 

— 

— 

16.76 + 0.51 

— 

— 

A 

— 

Sourcel 

03:27:52.002 

-29:09:09.65 

— 

— 

— 

14.55 + 0.22 

— 

— 

c 

— 

Source2 

03:27:52.025 

-29:09:12.14 

— 

— 

— 

7.47 + 0.14 

— 

— 

X 

J033210.8-270535 

HECDFS04 

03:32:10.840 

-27:05:34.18 

56 + 6 

61 + 6 

55 + 6 

14.57 + 0.26 

0.15 

0.44 X 0.35 

— 

— 

SourceO 

03:32:10.905 

-27:05:32.87 

— 

— 

— 

11.91 + 0.24 

— 

— 

c 

— 

Sourcel 

03:32:10.729 

-27:05:36.22 

— 

— 

— 

2.66 + 0.11 

— 

— 

c 

J033317.9-280907 

HECDFS13 

03:33:18.017 

-28:09:07.52 

95 + 6 

89 + 6 

63 + 6 

15.36 + 0.27 

0.14 

0.44 X 0.35 

— 

— 

SourceO 

03:33:18.006 

-28:09:07.55 

— 

— 

— 

10.11 + 1.30 

— 

— 

X 

— 

Sourcel 

03:33:18.032 

-28:09:07.39 

— 

— 

— 

5.25 + 1.37 

— 

— 

X 

J043340.5-540337 

HADFS04 

04:33:40.450 

-54:03:39.51 

74 + 6 

93 + 6 

84 + 6 

18.12 + 0.44 

0.19 

0.54 X 0.46 

— 

— 

SourceO 

04:33:40.455 

-54:03:40.29 

— 

— 

— 

9.25 + 0.30 

— 

— 

c 

— 

Sourcel 

04:33:40.501 

-54:03:40.05 

— 

— 

— 

6.09 + 0.33 

— 

— 

c 

— 

Source2 

04:33:40.472 

-54:03:38.33 

— 

— 

— 

2.78 + 0.19 

— 

— 

c 

J043619.3-552425 

HADFS02 

04:36:19.702 

-55:24:25.01 

102 + 6 

97 + 6 

81 + 5 

16.79 + 0.40 

0.19 

0.54 X 0.46 

— 

— 

SourceO 

04:36:19.706 

-55:24:24.41 

— 

— 

— 

7.81 + 0.47 

— 

— 

X 

— 

Sourcel 

04:36:19.698 

-55:24:25.27 

— 

— 

— 

8.99 + 0.58 

— 

— 

X 

J043829.7-541831 

HADFSll 

04:38:30.883 

-54:18:29.38 

19 + 6 

39 + 5 

52 + 6 

28.47 + 0.64 

0.19 

0.54 X 0.46 

— 

— 

SourceO 

04:38:30.780 

-54:18:31.79 

— 

— 

— 

21.19 + 0.51 

— 

— 

c 

— 

Sourcel 

04:38:30.970 

-54:18:26.60 

— 

— 

— 

7.28 + 0.30 

— 

— 

c 

J044103.8-531240 

HADFSIO 

04:41:03.942 

-53:12:41.01 

47 + 6 

58 + 6 

58 + 6 

17.44 + 0.39 

0.20 

0.55 X 0.45 

— 

— 

SourceO 

04:41:03.866 

-53:12:41.33 

— 

— 

— 

9.61 + 0.25 

— 

— 

X 

— 

Sourcel 

04:41:04.000 

-53:12:40.10 

— 

— 

— 

4.59 + 0.23 

— 

— 

X 

— 

Source2 

04:41:03.912 

-53:12:42.09 

— 

— 

— 

3.24 + 0.19 

— 

— 

X 

J044153.9-540350 

HADFSOl 

04:41:53.880 

-54:03:53.48 

76 + 6 

100 + 6 

94 + 6 

32.79 + 0.47 

0.19 

0.54 X 0.45 

A 

J044946.9-525424 

HADFS09 

04:49:46.448 

-52:54:26.95 

98 + 6 

102 + 6 

72 + 6 

15.52 + 0.59 

0.19 

0.54 X 0.45 

— 

— 

SourceO 

04:49:46.603 

-52:54:23.66 

— 

— 

— 

8.24 + 0.26 

— 

— 

X 

— 

Sourcel 

04:49:46.301 

-52:54:30.26 

— 

— 

— 

4.86 + 0.34 

— 

— 

X 

— 

Source2 

04:49:46.280 

-52:54:26.06 

— 

— 

— 

2.42 + 0.35 

— 

— 

X 

J045026.5-524127 

HADFS08 

04:50:27.453 

-52:41:25.41 

142 + 6 

133 + 6 

90 + 6 

14.18 + 0.50 

0.19 

0.54 X 0.45 

— 

— 

SourceO 

04:50:27.092 

-52:41:25.62 

— 

— 

— 

6.17 + 0.28 

— 

— 

c 

— 

Sourcel 

04:50:27.806 

-52:41:25.10 

— 

— 

— 

8.01 + 0.43 

— 

— 

c 

J045057.5-531654 

HADFS03 

04:50:57.715 

-53:16:54.42 

119 + 6 

102 + 6 

63 + 6 

11.50 + 0.39 

0.19 

0.54 X 0.45 

— 

— 

SourceO 

04:50:57.610 

-53:16:55.09 

— 

— 

— 

7.12 + 0.22 

— 

— 

c 

— 

Sourcel 

04:50:57.805 

-53:16:56.96 

— 

— 

— 

2.12 + 0.14 

— 

— 

c 

— 

Source2 

04:50:57.741 

-53:16:54.54 

— 

— 

— 

2.27 + 0.26 

— 

— 

c 

J100056.6+022014 

HCOSMOS02 10:00:57.180 

+02:20:12.70 

70 + 6 

85 + 6 

71 + 6 

14.61 + 0.66 

0.15 

0.63 X 0.50 

— 

— 

SourceO 

10:00:56.946 

+02:20:17.35 

— 

— 

— 

5.26 + 0.26 

— 

— 

X 

— 

Sourcel 

10:00:57.565 

+02:20:11.26 

— 

— 

— 

3.77 + 0.32 

— 

— 

X 

— 

Source2 

10:00:56.855 

+02:20:08.93 

— 

— 

— 

1.69 + 0.25 

— 

— 

X 
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Table 1 — Continued 


lAU address^ Short name 

RAgyo 

(J2000) 

Dec870 

(J2000) 

*8250 

(mjy) 

‘8350 

(mJy) 

*8500 

(mJy) 

Ssro CTALMA 

(mJy) (mJy) 

^ALMA 
(" X ") 

Lens 

grade^ 

— Sources 

10:00:57.274 

+ 02 : 20 : 12.66 

— 

— 

— 

1.66 + 0.21 — 

— 

X 

— Sourced 

10:00:57.400 

+02:20:10.83 

— 

— 

— 

2.23 + 0.41 — 

— 

X 

JI00144.I+0257I2 HCOSMOSOI 

10:01:44.182 

+02:57:12.47 

86 + 6 

96 + 6 

71 + 6 

15.35 + 0.25 0.14 

0.64 X 0.49 

A 


^ lAU name = IHerMES S250 + lAU address 

^ A = strongly lensed, C = weakly lensed, X = unlensed. Discussion 
of lens grades are given in Section 3.2. 

3. MODEL FITS 

3.1. Model Fitting Methodology 

An interferometer measures visibilities at discrete 
points in the uv plane. This is why pixel-to-pixel er¬ 
rors in the inverted and deconvolved surface brightness 
map of an astronomical source are correlated. The best 
way to deal with this situation is to compare model and 
data visibilities rather than surface brightness maps. The 
methodology used in this paper is similar in many as¬ 
pects to that used in Bussmann et al. (2012), who pre¬ 
sented the first lens model derived from a visibility-plane 
analysis of interferometric imaging of a strongly lensed 
DSFG discovered in wide-field submm surveys as well as 
in Bussmann et al. (2013), where this work was extended 
to a statistically significant sample of 30 objects. It also 
bears some resemblence to the method used in Hezaveh 
et al. (2013), who undertake lens modeling of interfero¬ 
metric data in the visibility plane. We summarize impor¬ 
tant information on the methodology here, taking care to 
highlight where any differences occur between this work 
and that of our previous efforts. 

We created and made publicly available custom soft¬ 
ware, called UVMCMCFIT, which is capable of modeling 
all of the ALMA sources in this paper efficiently. 

Sources are assumed to be elliptical Gaussians that are 
parameterized by the following six free parameters: the 
position of the source (relative to the primary lens if a 
lens is present, Ao^s and A^), the total intrinsic flux 
density the effective radius length {r^ = ^/o^), 

the axial ratio = h^/a^)^ and the position angle (^g? 
degrees east of north). The use of an elliptical Gaussian 
represents a simplification from the Sersic profile (Sersic 
1968) that is permitted based on the relatively weak con¬ 
straints on the Sersic index found in our previous work 
(Bussmann et al. 2012, 2013). 

When an intervening galaxy (or group of galaxies) is 
present along the line of sight, UVMCMCFIT accounts for 
the deflection of light caused by this structure using a 
simple ray-tracing routine that is adopted from a Python 
routine written by A. Bolton This represents a signif¬ 
icant difference from Bussmann et al. (2012) and Buss¬ 
mann et al. (2013), where we used the publicly available 
Gravlfns software (Keeton 2001) to map emission from 
the source plane to the image plane for a given lens- 
ing mass distribution. Gravlfns has a wide range of 
lens mass profiles as well as a sophisticated algorithm for 
mapping source plane emission to the image plane, but it 
also comes with a significant input/output penalty that 

http://www.physics.utah.edu/~bolton/python_lens_demo/ 


makes parallel computing prohibitively expensive. For 
example, modeling a simple system comprising one lens 
and one source typically required 24-48 hours using the 
old software, whereas the same system can be modeled 
in less than one hour with the pure-Python code (tests 
of the Bolton ray-tracing routine indicate it produces 
results consistent with Gravlfns). The use of pure- 
Python code for tracing the deflection of light rays is a 
critical component of making UVMCMCFIT computation¬ 
ally feasible. 

In UVMCMCFIT, lens mass profiles are represented by 
Aliens singular isothermal ellipsoid (SIE) profiles, where 
Mens is the number of lensing galaxies found from the 
best available optical or near-IR imaging (a multitude of 
evidence supports the SIE as a reasonable choice; for a 
review, see Treu 2010). Each SIE is fully described by 
the following five free parameters: the position of the 
lens on the sky relative to the arbitrarily chosen “image 
center” based on the ALMA 870 jivn emission and any 
lensing galaxies seen in the optical or near-IR (Aliens 
and Allens; these can be compared with the position of 
the optical or near-IR counterpart relative to the “image 
center”: Aq^nir and A^nir), the mass of the lens (pa¬ 
rameterized in terms of the angular Einstein radius, ^e), 
the axial ratio of the lens (giens = ^lens/niens), and the 
position angle of the lens (0iens; degrees east of north). 
Unless otherwise stated, when optical or near-IR imag¬ 
ing suggests the presence of additional lenses (see Eig- 
ure 3), we estimate centroids for each lens by eye and 
fix the positions of the additional lenses with respect to 
the primary lens. Each additional lens thus has three 
parameters: 6 >e, giens, and ^Rns- 

The total number of free parameters for any given sys¬ 
tem is — 5 T 3 X (^N\ens 1) T 6 * Agource? whcrc 

Asource IS the uumbcr of Gaussian profiles used. 

We assume secondary, tertiary, etc., lenses are located 
at the same redshift as the primary lens. If this assump¬ 
tion were incorrect, to first order only the conversion 
from an angular Einstein radius to a physical mass of 
the lensing galaxy would be affected. As the physical 
masses of the lensing galaxies are not the focus of this 
work, this assumption is reasonable. 

We use uniform priors for all model parameters. The 
prior on the position of the lenses covers ±076 (170) 
in both RA and Dec, a value that reflects the 1-a 
absolute astrometric solution between the ALMA and 
optical/near-IR images of 0'.'2 (074) for SDSS-based 
(2MASS-based) astrometric calibration. In Section 3.2, 
we discuss the level of agreement between the astrom¬ 
etry from the images and the astrometry from the lens 
modeling on an object-by-object basis. Eor 6 >e, the prior 
covers 07l — 6". The axial ratios of the lenses and sources 
are restricted to be giens > 0.3 and > 0.2. This as¬ 
sumption is justified for the lenses because our optical 
observations reveal lenses that are not highly elliptical 
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31.50s 31.00s 30.50s 2h18m30.00s 



02.00s 01.50s 2h22m01.00s 






Figure 2. ALMA 870/xm images (color scale, units of mJybeam”^) of HerMES DSFGs (images have not been corrected for primary 
beam attenuation). Contours (black and white) trace 250p,m or 500p,m emission from Herschel (starting at 4(7 and increasing by factors 
of 2, where a = 7 mJy). North is up, east is left. The FWHM size of the ALMA synthesized beam is shown in the lower left corner of each 
panel. A solid white circle shows the FWHM size of the primary beam. Dashed squares identify the regions of each image that are shown 
in greater detail in Fig. 3. 
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Figure 2. Continued. 
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Figure 2. Continued. 

and because we expect dark matter to be more spheri¬ 
cally distributed than the stars in leasing galaxies. For 
the sources, our ellipticity limit is primarily designed to 
aid numerical stability in the lens modeling. No prior is 
placed on the position angle of the lens or source. The 
intrinsic flux density for any source is allowed to vary 
from 0.1 mJy to the total flux density observed by ALMA 
(we ensure that the posterior PDF of the intrinsic flux 
density shows no signs of preferring a value lower than 
0.1 mJy). The source position is allowed to vary over any 
reasonable range necessary to fit the data (typically, this 
is ±l"-2"). The effective radius is allowed to vary from 
0'.'01-T.'5. 

The surface brightness map generated as part of UVM- 
CMCFIT is then converted to a “simulated visibility” 
dataset (Idnodei) in much the same way as MIRIAD’s 
UVMODEL routine. Indeed, the code used in UVMCM- 
CFIT is a direct Python port of UVMODEL (the use of 
UVMODEL itself is not possible for the same reason as 
Gravlens: constant input/output makes parallel com¬ 
puting prohibitively expensive), uvmcmcfit computes 
the Fourier transform of the surface brightness map and 
samples the resulting visibilities in a way that closely 
matches the sampling of the actual observed ALMA vis¬ 
ibility dataset (Valma)- 

The goodness of fit for a given set of model parameters 
is determined from the maximum likelihood estimate L 
according to: 

^ ^ ^ ( UALMA-Lnodell^ ^ log(27ra2)) (1) 

where a is the Icr uncertainty level for each visibility and 
is determined from the scatter in the visibilities within 


a single spectral window (this is a natural weighting 
scheme). 

We use EMCEE (Foreman-Mackey et al. 2013) to sam¬ 
ple the posterior probability density function (PDF) of 
our model parameters, emcee is a Markov chain Monte 
Carlo (MCMC) code that uses an affine-invariant ensem¬ 
ble sampler to obtain significant performance advantages 
over standard MCMC sampling methods (Goodman & 
Weare 2010). 

We employ a “burn-in” phase with 512 walkers and 
500-1000 iterations (i.e., ^ 250, 000 — 500, 000 samplings 
of the posterior PDF) to identify the best-fit model pa¬ 
rameters. This position then serves as the basis to initial¬ 
ize the “final” phase with 512 walkers and 10 iterations 
(i.e., 5,120 samplings of the posterior PDF) to determine 
uncertainties on the best-fit model parameters. 

During each MCMC iteration, we also measure the 
magnification factor at 870/im, /igTO? for each source. 
This is done simply by taking the ratio of the total 
flux density in the lensed image of the model (5'out) to 
the total flux density in the unlensed, intrinsic source 
model {Sin)- The use of an aperture when computing 
fiS 7 o is important when source profiles are used with sig¬ 
nificant flux at large radii (e.g., some types of Sersic 
profiles). For an elliptical Gaussian, such a step is un- 
neccessary (note that we did test this and found only 
^ 10% difference between figjo computed with and with¬ 
out an aperture). The best-fit value and la uncertainty 
on fis7o are drawn from the posterior PDF, as with the 
other parameters of the model. Exceptions are made 
for cases of weakly lensed sources where we have only 
upper limits on the Einstein radius (and hence upper 
limits on /igTo)- Ir such instances, we re-compute /igTO 
as the arithmetic mean of the limiting /igTO and unity: 
M870 = (/RSTO-Iimit + l)/2. The uncertainty in /igTO is 
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assumed to be equal to /i 870 -iimit ~ l)/2- 

Finally there are some important caveats to our ap¬ 
proach. The spatial resolution of the ALMA observa¬ 
tions is ^ 0''45, which is nearly always sufficient to re¬ 
solve the images of the lensed galaxy, but not always 
sufficient to resolve the images themselves. For this rea¬ 
son, in some cases the lens models may moderately over¬ 
predict the intrinsic sizes of the lensed galaxies and hence 
under-predict the magnification factors. In addition, our 
Gemini-S optical imaging may have missed optically faint 
lenses due to being at high redshift or dust-obscured (but 
not sufficiently active to be detected by ALMA). 

3.2. Individual Model Fits 

In this section, we present our model fits (as shown in 
Figure 3) and describe each source in detail. 

HELAISS02: Four sources are detected by ALMA, 
all of which are weakly lensed by a foreground galaxy 
seen in the HST image. To estimate the maximal magni¬ 
fication factors, we assume an Einstein radius of 1'.'5 for 
the lens (larger values predict counter images that are 
not seen by ALMA). The ALMA sources are all detected 
by IRAC and their mid-IR colors are similar, suggesting 
that they lie at the same redshift (see Figure 4). 

HXMM02: One source is detected by ALMA, and it 
is strongly lensed by one foreground galaxy seen in the 
HST image. The lensed source is not detected in the 
HST image. This object was first detected by Ikarashi 
et al. (2011) and also has high quality SMA imaging and 
an accompanying lens model that produces consistent 
results with those given here (Bussmann et al. 2013). 

HXMM31: Two sources are detected by ALMA, nei¬ 
ther of which are lensed. The faint, diffuse emission seen 
in the CFHT i-band image is atypical of leasing galax¬ 
ies. The nearest bright galaxy seen at i-band is located 
~ 18" southeast of the ALMA sources. 

HXMM29: Two sources are detected by ALMA, none 
of which appear to be lensed. The brighter ALMA source 
is weakly detected in the CFHT i-band image. 

HXMM07 : One source is detected by ALMA, and it 
is strongly lensed by one foreground galaxy detected in 
the Gemini-S image. There is a ~ 0"5 offset in the posi¬ 
tion of the foreground galaxy between the lens model and 
the Gemini-S image. Given the absolute astrometric rms 
uncertainty of 0"2 (based on SDSS), we do not consider 
this offset to be significant. The presence of a handful 
of ± 3(7 peaks in the residual map is likely an indication 
that our assumption of a single Gaussian to describe the 
source morphology is an oversimplification. 

HXMM20: Five sources are detected by ALMA, none 
of which appear to be lensed. There are a few faint 
smudges seen in the HST image which are likely to be 
the rest-frame optical counterparts to the ALMA sources. 
The ALMA sources are all arranged in a chain like shape, 
possibly suggestive of a larger filamentary overdensity in 
which they might reside. IRAC imaging provides support 
for this hypothesis (see Figure 4), as all of the ALMA 
sources are detected and have similar mid-IR colors. 

HXMMOl: Three sources are detected by ALMA, all 
of which are weakly lensed by two foreground galaxies 
seen in the HST and Keck/NIRC2 imaging. The ALMA 
imaging is broadly consistent with SMA data originally 
presented in Fu et al. (2013), with two bright sources 
and a much fainter third source very close to the more 


southern bright source. We assume Einstein radii of 0"5 
for both lenses in order to reproduce the approach used 
in Eu et al. (2013). This results in magnification factors 
for the three sources of /igTo ~ 1.6 — 1.7, consistent with 
Eu et al. (2013). 

HXMM04: One source is detected by ALMA, and it 
is weakly lensed by a foreground galaxy seen in the HST 
image. We assume an Einstein radius of 075 to represent 
the lensing scenario with maximum amplification. Due to 
the elliptical nature of the lens, this results in a maximum 
magnification factor of /igTO-iimit = 3.72±0.42. The HST 
morphology is complex: diffuse emission to the north of 
the lens could be a detection of the background source 
or could be a long spiral arm associated with the lensing 
galaxy. 

HXMM09: Two sources are detected by ALMA, both 
of which are weakly lensed by a single foreground galaxy 
detected in the HST image. An Einstein radius of 175 
is used to represent the “maximal lensing” scenario and 
results in maximal magnification factors of /igTO-Umit = 
2.25 ± 0.17 and yUgTO—limit = 1.48 ± 0.09. 

HXMM03: Three sources are detected by ALMA, all 
of which are weakly lensed by a foreground galaxy de¬ 
tected in the HST image and located ~ 6" from the 
ALMA sources. The central source is much brighter 
than the other two sources, which makes fitting a model 
challenging. We forced the positions of the second and 
third sources to be at least 075 and —075 away from 
the first source in declination, respectively. Eurther- 
more, we fixed the position of the lens to be located 
275 west and 075 south of the image centroid given in 
Table 1. We also fixed the Einstein radius to be l70, 
a typical value for isolated galaxies in this sample and 
in Bussmann et al. (2013). Because the source is so far 
from the lens, the maximal magnification factor is only 

MSTO —limit — 2.0 ± 0.1. 

HXMMll: Two sources are detected by ALMA, both 
of which are weakly lensed. This system is similar to 
HADES08, although the two ALMA sources are much 
closer and the lens must be less massive in order to 
avoid producing multiple images of the closest ALMA 
source. The fainter ALMA source has a much lower 
maximal magnification factor than the brighter source 
(M870-limit = 1.10 ± 0.01 VS. /i870-limit = 1.63 ± 0.11). 
Both ALMA sources are detected by IRAG and have sim¬ 
ilar mid-IR colors, suggesting they lie at similar redshifts 
(see Eigure 4). 

HXMM23: One source is detected by ALMA, and it 
is coincident (within the astrometric uncertainty) with 
a late-type galaxy seen in the HST image. Here, we 
assume that the HST source is the true counterpart to 
the ALMA source, implying that no lensing is occuring. 
Gonsistent with this hypothesis is that the SPIRE pho¬ 
tometry show blue colors that suggest this object is at 
low redshift. Note that models in which the late-type 
galaxy is lensing the ALMA source by a modest amount 
(M 870 <1.2) cannot be ruled out with the present data. 

HXMM22: One source is detected by ALMA, and it 
appears to be unlensed. A faint smudge seen in the HST 
image of this source is due to a star located 375 northeast 
of the ALMA source. 

HXMM05: One source is detected by ALMA, and 
it is weakly lensed by two foreground galaxies seen in 
the HST images. To compute the maximum magnifica- 
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tion factor, we assume an Einstein radius of 1" for the 
foreground lenses and fix the positions of both lenses ac¬ 
cording to the location of the foreground galaxies in the 
HST image. 

HXMM30: One source is detected by ALMA, and it 
is strongly lensed by one foreground galaxy detected in 
the Gemini-S image. As with HXMM07, there is a ~ 
0'.'5 offset between the lens position according to the lens 
model and the Gemini-S image. We do not consider this 
offset significant. An alternative model in which the lens 
is sub-mm luminous cannot be ruled out, but we consider 
this unlikely for a number of reasons. First, it is a more 
complex model (having two sources and one lens, rather 
than one source and one lens). Second, lenses are very 
rarely detected in sub-mm imaging. Third, the shape and 
location of the ALMA sources relative to the Gemini-S 
source are typical of strongly lensed objects (consistent 
with the very low residuals). Fourth, the alternative lens 
model predicts the lensed source to have an intrinsic flux 
density of ^ 13 mJy, which would make it the brightest 
source in the sample. 

HXMM12: One source is detected by ALMA, and it 
is weakly lensed by a group of foreground galaxies seen 
in the HST image. We assume an Einstein radius of 0'.'2 
for the nearest leasing galaxy and allow a ±0'.'4 (i.e., 2a) 
shift in its position relative to that indicated by the HST 
image (which has its astrometry tied to SDSS). We repre¬ 
sent the remaining members of the group as a single SIS 
(assumed to be spherical to simplify the model) located 
4'.'5 south and 4'.'5 east of the image centroid and having 
an Einstein radius of 2'.'0. This SIS is justified by the 
presence of several sources in this region of the HST im¬ 
age (not shown in Figure 3). This is meant to represent 
the “maximal leasing” scenario. The presence of two 3cr 
peaks located near the center of the residual image indi¬ 
cates that the model does not fit the data perfectly. This 
could be an indication that either of our assumptions for 
the lens potential or source structure are oversimplifica¬ 
tions. Higher resolution imaging is needed to determine 
the most likely cause. 

HECDFS12: This is a complex, well-constrained sys¬ 
tem. Two background sources are detected by ALMA: 
one is multiply imaged and the other is singly imaged. 
In addition, the lens is detected by ALMA (this is one of 
two sources in the entire Herschel-ALMA sample that is 
unresolved by ALMA). These facts work together to pro¬ 
vide very tight constraints on the system. Since the lens 
is detected by ALMA, its position relative to the lensed 
images is unambiguous. Also, because there is a strongly 
lensed source with multiple images, the Einstein radius 
of the lens is unambiguous. Finally, this source is de¬ 
tected (and unresolved) in the NRAO VLA Sky Survey 
(Gondon et al. 1998), having 5 'i.4ghz = 21.8 ± 0.8 mJy. 
Assuming all of this radio emission originates from the 
lens, this implies a spectral slope of a = —0.24 and is 
consistent with non-thermal emission from the lens. For 
this target, we show VIDEO imaging (Jarvis et al. 
2013). 

HECDFS04: Two sources are detected by ALMA, 
both of which are weakly lensed by a foreground galaxy 
seen in the HST image. There is also a 3cr peak coin¬ 
cident with an HST source that may be an indication 
that the lens has been detected by ALMA. We do not 
attempt to model this 3cr peak. We assume an Einstein 


radius of 0'.'5 for the lens, since larger values predict the 
existence of counter images that are not seen by ALMA. 
The second ALMA source is located ~ 5" from the lens 
and experiences a small but significant magnification of 
M870-iimit = 1.12 ± 0.02. Both ALMA sources appear 
to be detected by IRAG and have similar mid-IR colors, 
suggesting they he at the same redshift (see Figure 4). 

HECDFS13: This system is similar to HADFS02 
(mentioned below), except that here the two ALMA 
sources are separated by ^ 0'.'4 rather than 0'.'8 and one 
source is brighter than the other by a factor of 2. As¬ 
suming the two sources have similar mass-to-light ratios, 
their brightness ratios indicate major merger rather than 
minor merger activity. The projected physical distance is 
^2 — 3 kpc, assuming a redshift of z = 2 for the ALMA 
sources. This could be an example of a major merger 
approaching final coalescence and experiencing a signifi¬ 
cant boost in star-formation due to enhancements in the 
local gas density brought about by tidal forces during the 
merger. 

HADFS04: Three sources are detected by ALMA, all 
of which are weakly lensed by a foreground galaxy seen 
in the HST image. We assume an Einstein radius of 075 
for the lens, as values larger than this produce multiple 
images of the ALMA sources. Values for the Einstein 
radius that are smaller than 0'.'5 are unlikely based on 
the brightness of the lens, so the results we report for 
this object should be robust. 

HADFS02: Two sources are detected by ALMA. The 
nearest possible lens is located ~ 8" from the ALMA 
sources, indicating that leasing is likely to be irrelevant in 
this system. The two ALMA sources are similarly bright 
('S'870 = 8-27 ± 0.53 mJy and Ssto = 9.07 ± 0.27 mJy) 
and separated by ^ 078, corresponding to a projected 
physical distance of ^ 6 kpc. This distance is typical of 
the pericentric passage distance in both hydrodynami- 
cal simulations of major mergers (e.g., Hayward et al. 
2012a) and observations of major mergers (e.g., Ivison 
et al. 2007; Tacconi et al. 2008; Engel et al. 2010; Riechers 
et al. 2011b; Ivison et al. 2011). Two plausible scenarios 
are that HADFS02 represents a major merger that just 
experienced a first pass or is approaching final coales¬ 
cence, either of which would significantly enhance star- 
formation in the system. 

HADFSll: Two sources are detected by ALMA, 
both of which are weakly lensed by a group of small 
galaxies detected in the HST image. To estimate the 
maximum magnification factor, we represent the gravi¬ 
tational potential of the group with a single SIE lens and 
an Einstein radius of I'.'O. Values larger than this pro¬ 
duce additional counter images that are not seen in the 
ALMA imaging. 

HADFSIO: Three sources are detected by ALMA. 
We assume that all three are unlensed. There is a group 
of three sources detected in our Gemini-S optical imaging 
located ~ 1" east of the ALMA sources. This distance 
is so large that plausible mass ranges for the Gemini-S 
sources would imply at most a factor of 1.1 —1.2 boost 
in the apparent flux densities of the ALMA sources. We 
also tested a single-lens, single-source model in which the 
source is triply-imaged in the manner that is observed. 
The lens in this hypothetical model has an Einstein ra¬ 
dius of 1'.'2, requiring a very high mass to light ratio or 
a very high lens redshift to be consistent with the non- 
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detection in the Gemini-S data. Deep near-IR imaging 
is needed to confirm that this target is unlensed. 

HADFSOl: This is a single source that is strongly 
lensed by a foreground galaxy seen in the HST image. 
The lensed source is not detected by HST. The source is 
highly elongated {q^ = 0.31 ±0.01), but fits the data very 
well. The position of the lens according to the lens model 
is consistent with the position in the HST image, given 
the 0'.'4 fundamental uncertainty due to using the 2MASS 
system as the fundamental basis for the astrometry. 

HADFS09: Three sources are detected by ALMA, 
none of which appear to be lensed (the closest bright HST 
source is located ~ 13" away from the ALMA sources). 

HADFS08: Two sources are detected by ALMA, 
both of which are weakly lensed by a foreground galaxy 
in the HST image. The ALMA sources have the largest 
separation of any in our sample overall, around 10". We 


assume an Einstein radius of 1"5 for the foreground lens 
as a “maximal lensing” scenario. This results in maxi¬ 
mum magnification factors of //870-iimit = 2.3 ±0.1 and 
M870-iimit = 1.2 ± 0.1 for the two sources. 

HADFS03: Three sources are detected by ALMA, 
each of which is weakly lensed by a single bright fore¬ 
ground galaxy seen in the HST image. Alternative sce¬ 
narios involving strong lensing can be ruled out by the 
location of the lens: ~ 2" — 3" north of the centroid of 
the ALMA sources (the rms error in the astrometry is set 
from 2MASS at a level of ~ 0"4) as well as the atypical 
location and fluxes of the ALMA sources relative to each 
other. To obtain the maximum magnification factor, we 
assume an Einstein radius of 0"5 and fix the position 
angle of the lens to be between 40°-50° to match the 
orientation seen in the HST image. Larger Einstein radii 
can be ruled out by the absence of counter images north 
of the lens. 


Table 2 

Lens properties from parameters of model fits to ALMA sources (parameters 
are described in Section 3.1). Parameters without uncertainties were fixed to 

the given value. 


Short name 

ARAsyo 

(") 

ADec87o 

(") 

(") 

^lens 

*^lens 

(deg) 

HELAISS02.Lens0 

-1.59 ±0.20 

2.25 ±0.19 

1.500 

0.790 ± 0.067 

44± 16 

HXMM02.Lens0 

0.01 ±0.01 

-0.24 ±0.01 

0.507 ± 0.004 

0.596 ± 0.009 

157 ± 10 

HXMM07.Lens0 

-0.27 ±0.03 

0.04 ±0.13 

0.928 ± 0.007 

0.902 ± 0.024 

26 ±7 

HXMMOl.LensO 

2.05 

0.60 

0.500 

0.801 ± 0.062 

48 ±14 

HXMMOl.Lensl 

-2.80 

1.00 

0.500 

0.882 ± 0.072 

90± 17 

HXMM04.Lens0 

0.17 ±0.03 

0.04 ± 0.03 

0.500 

0.547 ± 0.050 

11 ± 16 

HXMMOg.LensO 

1.40 ±0.07 

0.19 ±0.05 

1.000 

0.663 ± 0.094 

64± 16 

HXMMOS.LensO 

-2.50 

-0.50 

1.000 

1.000 

0 

HXMMll.LensO 

0.82 ±0.12 

2.95 ±0.10 

0.500 

0.706 ±0.124 

67± 11 

HXMM05.Lens0 

2.80 

-1.40 

1.000 

0.531 ±0.180 

45 ±14 

HXMMOS.Lensl 

-1.90 

2.50 

1.000 

0.569 ±0.197 

67± 16 

HXMM30.Lens0 

-0.03 ± 0.02 

0.05 ±0.01 

0.743 ± 0.008 

0.703 ± 0.050 

26 ±10 

HXMM12.LensO 

-0.22 ±0.20 

-0.25 ± 0.24 

0.200 

0.672 ± 0.090 

30± 16 

HXMM12.Lensl 

4.50 

-4.50 

2.000 

1.000 

0 

HECDFS12.LensO 

0.22 

-1.75 

1.354 ±0.006 

0.955 ± 0.007 

80± 16 

HECDFS04.Lens0 

1.01 ±0.02 

2.10 ±0.01 

0.500 

0.807 ± 0.006 

176 ± 13 

HADFS04.Lens0 

-0.56 ±0.13 

0.11 ±0.07 

0.500 

0.662 ±0.135 

37± 12 

HADFSll.LensO 

0.41 ±0.04 

0.27 ±0.12 

1.000 

0.723 ± 0.068 

82 ±19 

HADFSOl.LensO 

-0.19 ±0.01 

0.25 ±0.01 

1.006 ±0.004 

0.794 ± 0.008 

99 ±10 

HADFSOS.LensO 

-3.59 ±0.06 

-2.32 ±0.06 

1.500 

0.897 ± 0.047 

74± 18 

HADFS03.Lens0 

-0.40 ±0.08 

1.32 ±0.06 

1.000 

0.707 ±0.141 

93 ±17 

HCOSMOSOl.LensO 

-0.12 ±0.01 

0.28 ±0.02 

0.956 ± 0.005 

0.775 ± 0.025 

72 ±10 


HCOSMOS02: Eive sources are detected by ALMA 
(the brightest of which was already known; Smolcic et al. 
2012 ), none of which appear to be lensed. Previous re¬ 
search has shown this to be an overdense region (this 
object is called C0SB03 in Aravena et al. 2010) with an 
optical and near-IR photometric redshift of z = 2.3 —2.4. 
Our ALMA imaging offers the first convincing evidence 
that the associated galaxies in the overdensity are sub- 
mm bright and thus intensely star-forming. There are a 
number of 2 — 3cr peaks in the map that could be real. 
This would further increase the multiplicity rate for this 
object, but we caution that there are also negative peaks 
of similar amplitude (i.e., 2 — 3cr) present in this map. 
Some of the ALMA sources have counterparts detected 
in the HST image, whereas all of the ALMA sources are 
detected by IRAC (see Eigure 4). Their mid-IR colors 
are similar, providing further evidence that the ALMA 


sources lie at the same redshift. 

HCOSMOSOl: This system is similar to HADESOl: 
a single source that is strongly lensed by a foreground 
galaxy seen in the HST image. In fact, the background 
source is also detected by HST as well as Keck/NIRC2 
adaptive optics imaging, and a lens model has been pub¬ 
lished based on these data (Calanog et al. 2014). The 
morphology of the lensed emission is very different be¬ 
tween the Keck and ALMA imaging, suggesting differen¬ 
tial magnification is important in this object. The very 
small sizes of the sources are consistent with this as well 
(rs = 0.023 ± 0.003", Keck and Ts = 0.055 ± 0.007", 
Alma). Adopting a redshift of z = 2 for the lensed 
source implies physical sizes of ~ 150 pc and ~ 300 pc 
for the rest-frame optical and rest-frame EIR emission, 
respectively. 
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Figure 3. Model fits for each target in the Herschel-ALMA sample, 3 panels per target. North is up and east is left. Left: ALMA 870p,m 
imaging (red contours, starting at ±3(7 and increasing by factors of 2) overlaid on best available optical or near-IR imaging (grayscale, 
with telescope and filter printed in upper right corner). The location and morphology of all sources used in the model are represented by 
magenta ellipses. If a lens is present, its location is given by a black circle and its critical curve is traced by an orange line. The FWHM 
size of the ALMA synthesized beam is shown in the lower left corner of each panel. Middle: Same as leff but showing best-fit model in 
grayscale. Numbers indicate the location of sources. Right: Same as lefh but showing residual image obtained from subtracting best-fit 
model from the data. 
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Figure 3. Continued. 
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Table 3 

Intrinsic properties from parameters of model fits to ALMA sources 
(parameters are described in Section 3.1). Uncertainties in flux densities do 
not include absolute calibration uncertainty of 10%. Note that some 
parameters such as source size may become unreliable when the signal to 
noise ratio is below 10 (e.g., Simpson et al. 2015a). 


Short name 

ARAsro 

(") 

ADecsro 

(") 

■S'870 

(mJy) 

rs 

(") 


qs 


(deg) 

M870 

HELAISS02.0 

3.113 ±0.160 

-3.112 ±0.155 

8.02 

± 

0.15 

0.096 

± 

0.005 

0.80 

± 

0.05 

91 ±6 

1.15 

± 

0.07 

HELAISS02.1 

-0.111 ±0.114 

-2.172 ±0.183 

3.24 

± 

0.12 

0.065 

± 

0.008 

0.84 

± 

0.05 

87 ±7 

1.34 

± 

0.17 

HELAISS02.2 

-1.470 ±0.158 

1.774 ±0.145 

3.27 

± 

0.25 

0.105 

± 

0.016 

0.86 

± 

0.04 

120 ±7 

1.27 

± 

0.13 

HELAISS02.3 

4.039 ±0.165 

-7.216 ±0.174 

2.22 

± 

0.18 

0.124 

± 

0.020 

0.79 

± 

0.05 

77 ±7 

1.08 

± 

0.04 

HXMM02.0 

-0.278 ± 0.008 

0.239 ±0.011 

11.85 

^± 

:0.11 

0.122 

± 

0.003 

0.64 

± 

0.02 

62 ±2 

5.33 

± 

0.19 

HXMM31.0 

-1.380 ±0.010 

-1.025 ±0.010 

6.79 

± 

0.37 

0.141 

± 

0.011 

0.80 

± 

0.12 

134 ± 36 


— 


HXMM31.1 

1.311 ±0.011 

1.124 ±0.010 

4.01 

± 

0.26 

0.070 

± 

0.018 

0.59 

± 

0.22 

52 ±56 


— 


HXMM29.0 

0.114 ±0.009 

0.451 ± 0.008 

5.46 

± 

0.30 

0.088 

± 

0.012 

0.82 

± 

0.14 

90 ±44 


— 


HXMM29.1 

-0.236 ±0.034 

-0.562 ± 0.030 

1.82 

± 

0.38 

0.116 

± 

0.051 

0.70 

± 

0.20 

88 ±55 


— 


HXMM07.0 

0.016 ±0.238 

-0.016 ±0.283 

3.43 

± 

0.07 

0.074 

± 

0.007 

0.32 

± 

0.02 

66 ±2 

8.49 

± 

1.13 

HXMM20.0 

-2.308 ±0.012 

-2.275 ±0.011 

7.15 

± 

0.44 

0.089 

± 

0.014 

0.63 

± 

0.16 

58 ±27 


— 


HXMM20.1 

0.828 ± 0.025 

-0.278 ± 0.023 

3.52 

± 

0.41 

0.137 

± 

0.026 

0.84 

± 

0.10 

74 ±44 


— 


HXMM20.2 

-0.211 ±0.017 

-1.647 ±0.014 

3.42 

± 

0.26 

0.058 

± 

0.020 

0.80 

± 

0.13 

84 ±45 


— 


HXMM20.3 

-1.505 ±0.157 

-1.981 ±0.064 

2.46 

± 

0.47 

0.283 

± 

0.198 

0.67 

± 

0.17 

81 ±21 


— 


HXMM20.4 

2.588 ±0.155 

2.611 ±0.218 

0.94 

± 

0.18 

0.459 

± 

0.246 

0.58 

± 

0.15 

96 ±51 


— 


HXMMOl.O 

-1.503 ±0.013 

0.395 ±0.017 

11.61 

. ± 

:0.23 

0.090 

± 

0.005 

0.56 

± 

0.06 

12 ±19 

1.39 

± 

0.19 

HXMMOl.l 

-2.563 ±0.018 

-1.337 ±0.017 

9.56 

± 

0.26 

0.116 

± 

0.006 

0.34 

± 

0.03 

2±1 

1.21 

± 

0.10 

HXMM01.2 

-2.622 ±0.025 

-0.552 ± 0.025 

1.45 

± 

0.20 

0.077 

± 

0.025 

0.66 

± 

0.18 

134 ± 33 

1.29 

± 

0.15 

HXMM04.0 

0.095 ±0.021 

0.442 ± 0.025 

8.49 

± 

0.20 

0.117 

± 

0.007 

0.52 

± 

0.07 

-2 ±5 

2.36 

± 

0.68 

HXMM09.0 

-0.392 ±0.039 

-0.740 ± 0.051 

5.51 

± 

0.19 

0.064 

± 

0.006 

0.42 

± 

0.06 

75 ±5 

1.62 

± 

0.31 

HXMM09.1 

-1.507 ±0.073 

0.805 ± 0.053 

5.14 

± 

0.15 

0.033 

± 

0.010 

0.46 

± 

0.18 

116 ± 14 

1.24 

± 

0.12 

HXMM03.0 

5.180 ±0.003 

0.924 ± 0.003 

12.25 

^± 

:0.24 

0.130 

± 

0.004 

0.53 

± 

0.03 

-25 ±2 

1.50 

± 

0.25 

HXMM03.1 

7.560 ± 0.023 

2.051 ±0.028 

1.46 

± 

0.13 

0.093 

± 

0.007 

0.73 

± 

0.11 

22 ±25 

1.50 

± 

0.25 

HXMM03.2 

7.663 ± 0.035 

0.755 ± 0.033 

1.35 

± 

0.12 

0.096 

± 

0.005 

0.73 

± 

0.11 

-11 ±37 

1.50 

± 

0.25 

HXMMll.O 

-0.844 ±0.111 

-0.648 ± 0.081 

6.24 

± 

0.24 

0.106 

± 

0.007 

0.26 

± 

0.03 

54 ±2 

1.31 

± 

0.16 

HXMMll.l 

-0.596 ±0.122 

-4.592 ± 0.098 

3.38 

± 

0.35 

0.168 

± 

0.023 

0.59 

± 

0.16 

139 ± 41 

1.05 

± 

0.03 

HXMM23.0 

0.101 ±0.011 

-0.050 ± 0.009 

2.93 

± 

0.15 

0.020 

± 

0.008 

0.68 

± 

0.20 

89 ±49 


— 


HXMM22.0 

-0.076 ± 0.004 

0.024 ± 0.004 

lO.lf 

)± 

:0.28 

0.085 

± 

0.010 

0.52 

± 

0.11 

152 ±6 


— 


HXMM05.0 

-3.505 ± 0.094 

1.937 ±0.081 

12.82 

;± 

:0.31 

0.095 

± 

0.006 

0.59 

± 

0.06 

142 ±5 

1.40 

± 

0.20 

HXMM30.0 

0.153 ±0.024 

-0.073 ±0.011 

0.84 

± 

0.01 

0.019 

± 

0.003 

0.20 

± 

0.00 

109 ±1 

27.1f 

► ± 

:4.61 

HXMM12.0 

1.520 ±0.168 

-0.683 ± 0.243 

9.91 

± 

0.24 

0.115 

± 

0.005 

0.72 

± 

0.07 

69 ±8 

1.57 

± 

0.29 

HECDFS12.0 

-0.348 ± 0.006 

0.077 ± 0.004 

2.02 

± 

0.06 

0.085 

± 

0.004 

0.38 

± 

0.03 

134 ±3 

8.29 

± 

0.19 

HECDFS12.1 

-0.342 ± 0.005 

2.489 ± 0.008 

11.54 ± 

:0.18 

0.147 

± 

0.003 

0.65 

± 

0.02 

14 ±2 

1.26 

± 

0.13 

HECDFS12.2 

0.000 

0.000 

7.47 

± 

0.14 

0.026 

± 

0.009 

0.79 

± 

0.15 

85 ±63 


— 


HECDFS04.0 

-0.011 ±0.011 

-0.347 ± 0.004 

6.02 

± 

0.12 

0.096 

± 

0.005 

0.35 

± 

0.03 

91 ±2 

1.98 

± 

0.49 

HECDFS04.1 

-2.366 ±0.024 

-3.752 ± 0.007 

2.51 

± 

0.10 

0.032 

± 

0.012 

0.68 

± 

0.19 

93 ±55 

1.06 

± 

0.03 

HECDFS13.0 

-0.156 ±0.011 

-0.034 ±0.011 

10.11 

.± 

: 1.30 

0.099 

± 

0.012 

0.52 

± 

0.12 

123 ±7 


— 


HECDFS13.1 

0.221 ±0.061 

0.127 ±0.018 

5.25 

± 

1.37 

0.109 

± 

0.024 

0.38 

± 

0.08 

88 ±7 


— 


HADFS04.0 

0.333 ±0.101 

-0.513 ± 0.040 

6.85 

± 

0.22 

0.091 

± 

0.006 

0.39 

± 

0.05 

142 ±4 

1.35 

± 

0.17 

HADFS04.1 

0.865 ±0.123 

-0.420 ± 0.041 

5.03 

± 

0.27 

0.165 

± 

0.013 

0.43 

± 

0.06 

141 ±4 

1.21 

± 

0.10 

HADFS04.2 

0.604 ±0.108 

0.739 ± 0.077 

1.99 

± 

0.14 

0.077 

± 

0.015 

0.75 

± 

0.16 

101 ± 40 

1.40 

± 

0.20 

HADFS02.0 

0.067 ±0.008 

0.588 ±0.015 

7.81 

± 

0.47 

0.136 

± 

0.012 

0.38 

± 

0.06 

23 ±5 


— 


HADFS02.1 

-0.060 ± 0.009 

-0.268 ±0.018 

8.99 

± 

0.58 

0.193 

± 

0.015 

0.42 

± 

0.06 

17 ±4 


— 


HADFSll.O 

-1.340 ±0.043 

-1.816 ±0.119 

17.51 

.± 

:0.42 

0.225 

± 

0.006 

0.46 

± 

0.02 

178 ±1 

1.21 

± 

0.11 

HADFSll.l 

0.658 ± 0.039 

1.569 ±0.111 

5.78 

± 

0.24 

0.180 

± 

0.010 

0.25 

± 

0.02 

167 ±2 

1.26 

± 

0.13 

HADFSIO.O 

-1.126 ±0.005 

-0.319 ± 0.004 

9.61 

± 

0.25 

0.073 

± 

0.010 

0.67 

± 

0.15 

133 ± 24 


— 


HADFSlO.l 

0.876 ±0.011 

0.908 ± 0.009 

4.59 

± 

0.23 

0.048 

± 

0.019 

0.71 

± 

0.19 

84 ±43 


— 


HADFS10.2 

-0.437 ±0.017 

-1.088 ±0.016 

3.24 

± 

0.19 

0.093 

± 

0.020 

0.58 

± 

0.20 

131 ± 38 


— 


HADFSOl.O 

0.131 ±0.005 

-0.105 ±0.006 

3.17 

± 

0.05 

0.128 

± 

0.005 

0.30 

± 

0.01 

24 ± 1 

10.34 ± 

:0.47 

HADFS09.0 

2.343 ± 0.007 

3.284 ± 0.005 

8.24 

± 

0.26 

0.109 

± 

0.008 

0.70 

± 

0.11 

92 ±14 


— 


HADFS09.1 

-2.191 ±0.013 

-3.320 ±0.011 

4.86 

± 

0.34 

0.099 

± 

0.019 

0.53 

± 

0.17 

135 ± 24 


— 


HADFS09.2 

-2.503 ±0.035 

0.886 ±0.019 

2.42 

± 

0.35 

0.122 

± 

0.040 

0.51 

± 

0.17 

89 ±20 


— 


HADFS08.0 

-0.868 ±0.050 

0.938 ± 0.048 

3.74 

± 

0.17 

0.055 

± 

0.010 

0.83 

± 

0.09 

131 ± 20 

1.65 

± 

0.32 

HADFS08.1 

7.496 ± 0.058 

2.190 ±0.059 

7.28 

± 

0.39 

0.179 

± 

0.012 

0.59 

± 

0.08 

63 ±6 

1.10 

± 

0.05 

HADFS03.0 

-0.734 ± 0.069 

-1.070 ±0.053 

5.39 

± 

0.17 

0.112 

± 

0.006 

0.41 

± 

0.05 

45 ±3 

1.32 

± 

0.16 

HADFS03.1 

1.415 ±0.056 

-2.912 ±0.059 

1.87 

± 

0.13 

0.059 

± 

0.018 

0.54 

± 

0.12 

94 ±48 

1.13 

± 

0.07 

HADFS03.2 

0.427 ±0.055 

-0.514 ±0.059 

1.22 

± 

0.14 

0.084 

± 

0.020 

0.50 

± 

0.13 

125 ± 14 

1.86 

± 

0.43 

HCOSMOS02.0 

-3.507 ±0.012 

4.659 ±0.013 

5.26 

± 

0.26 

0.073 

± 

0.017 

0.70 

± 

0.12 

94 ±34 


— 


HCOSMOS02.1 

5.780 ±0.019 

-1.434 ±0.026 

3.77 

± 

0.32 

0.094 

± 

0.029 

0.76 

± 

0.13 

106 ± 65 


— 


HCOSMOS02.2 

-4.869 ± 0.049 

-3.769 ± 0.050 

1.69 

± 

0.25 

0.198 

± 

0.051 

0.65 

± 

0.13 

72 ±41 


— 


HCOSMOS02.3 

1.410 ±0.031 

-0.035 ± 0.033 

1.66 

± 

0.21 

0.101 

± 

0.042 

0.71 

± 

0.13 

74 ±42 


— 


HCOSMOS02.4 

3.301 ±0.083 

-1.864 ±0.060 

2.23 

± 

0.41 

0.312 

± 

0.060 

0.67 

± 

0.13 

78 ±32 


— 


HCOSMOSOl.O 

0.136 ±0.011 

-0.220 ±0.016 

1.03 

± 

0.02 

0.068 

± 

0.006 

0.27 

± 

0.04 

164 ±2 

14.8e 

)± 

: 1.90 
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4. RESULTS 


4.1. De-lensing the ALMA Sample 

The combination of our optical or near-IR imaging and 
our deep, high-resolution ALMA imaging permits us to 
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Figure 4. ALMA 870 iim imaging (white contours, starting at Act and increasing by factors of 2) overlaid on color composite IRAC 
imaging (blue = 3.6p,m, green = 4.5/rm, red = 8.0/rm). All panels are 27" on a side. North is up and east is left. The synthesized beam 
is represented in the lower right corner of each panel. Each of the ALMA counterparts are detected in the IRAC imaging. In addition, 
the IRAC colors of ALMA sources are broadly consistent, providing some evidence that they are at the same redshift and not physically 
unassociated blends along the line of sight. 


take the first step towards mapping the foreground struc¬ 
ture along the line of sight to the ALMA sources. With 
such maps in hand for all of our targets, we can estimate 
the impact that leasing has on the intrinsic properties of 
the ALMA sources. In other words, we can “de-lens” the 
Herschel-PsljyiK sample. 

Figure 5 shows the observed (i.e., apparent) and intrin¬ 
sic (i.e., de-lensed) distributions of S'sro- Lensing has the 
strongest effect on S'gTO- fhe median fiux density in the 
Herschel-ALyiA sample drops by a factor of 1.6 when 
lensing is taken into account. A two-sided Kolmogorov- 
Smirnov (KS) test yields a p-value of 0.044, suggesting 
the apparent and intrinsic fiux density distributions are 
inconsistent with being drawn from the same parent pop¬ 
ulation. Even if strongly lensed sources are removed from 
the sample, the median intrinsic fiux density is 1.3 times 
lower than the median apparent fiux density. Remov¬ 
ing the unlensed sources from consideration pushes this 
factor back to 1.6. At these levels, failing to correct for 
amplification due to gravitational lensing will be a signif¬ 
icant source of error, since the absolute calibration un¬ 
certainty is typically of order 5-10%. When discussing 
the intrinsic properties of bright sources (including their 
number counts, e.g. Wyithe et al. 2011) discovered in 
wide-field FIR or mm surveys, it is critical to consider 
the effects of lensing. 

For comparison, we also show the cumulative distribu¬ 
tion of S'sTO for the ALESS sample (including the com¬ 
pleteness limit of LESS of 4.5 mJy). ALESS is the only 
existing sample of DSFGs with interferometric follow-up 
of a sensitivity and angular resolution that is compara¬ 
ble to our ALMA data, so it is the best sample with 
which to compare our results. The significant overlap 
in S'sTO between our sample and ALESS is evidence that 
the DSFGs in our sample have higher S'soo/'S'sto ratios 
(even when the effect of lensing in our sample is taken 
into account) than the DSFGs in ALESS (recall Figure 1, 
which shows that ALESS sources have much lower ^Ssoo 
than our targets). This difference is likely due to differ¬ 
ences in dust temperature and/or redshift distributions 
of the two samples and probably arises from selection 
effects. 

The effect on the other source parameters (rg, angu¬ 
lar separation (the angular distance between an ALMA 
source and the centroid of all the ALMA sources for a 
given Herschel DSFG), and q^) is less pronounced. The 
median source size decreases by a factor of 1.2 in the 
Herschel-ALMA sample after accounting for lensing, but 
the two-sided KS test reveals a p-value of 0.174, suggest¬ 


ing that we cannot rule out the null hypothesis that both 
size distributions were drawn from the same parent dis¬ 
tribution. We find no significant difference between the 
axial ratios of the apparent and intrinsic distributions, 
as well as between the angular separations of apparent 
and intrinsic distributions (two-sided KS test p-values of 
0.984 and 0.920, respectively). 

Finally, the brightest source in the Herschel-ALMA 
sample is HADES 11.0, with an intrinsic fiux density of 
Ss7o = 17.5 ± 0.4 mJy. However, there are also two ob¬ 
jects with multiple sources that have separations smaller 
than 1", which have summed fiux densities comparable 
to this; namely HADFS02 (16.8 mJy) and HEGDFS13 
(15.3 mJy). This is approaching the values found in the 
most extreme systems, such as GN20 (20.6 mJy, Pope 
et al. 2006) and HFLS3 (15-20 mJy; Riechers et al. 2013; 
Gooray et al. 2014; Robson et al. 2014). It is a level 
that is extremely difficult to reproduce in simulations 
(e.g., Narayanan et al. 2010). One possibility is that the 
objects with multiple sources represent blends of physi¬ 
cally unassociated systems. We explore this possibility 
via comparison to theoretical models in Section 4.3, but 
a direct empirical test requires redshift determinations 
for each source and is beyond the scope of this paper. 

4.2. Multiplicity in the ALMA Sample 

The second key result from our deep, high-resolution 
ALMA imaging is a firm measurement of the rate of mul¬ 
tiplicity in Herschel DSFGs. We find that 20/29 Herschel 
DSFGs break down into multiple ALMA sources, imply¬ 
ing a multiplicity rate of 69%. However, 5/9 of the single¬ 
component systems are strongly lensed. If these five are 
not considered, then the multiplicity rate increases to 
80%. Such a high rate of multiplicity is consistent with 
theoretical models (e.g., Hayward et al. 2013a, hereafter 
HB13). 

In comparison, the 69 DSFGs in the MAIN ALESS cat¬ 
alog show a multiplicity rate of 35 - 40% (Hodge et al. 
2013). Smoothing our ALMA images and adding noise 
to match the resolution and sensitivity of ALESS re¬ 
sults in a multiplicity rate of 55% (four objects with 
sources that are separated by < 1" become single sys¬ 
tems). On the other hand, the ALESS sources are 
much fainter overall, having a median 870 /im fiux den¬ 
sity of Ss 7 o ~ 6mJy, compared to S'gTO = 14.9 mJy 
in our Herschel-ALMA sample. Thus, the evidence fa¬ 
vors brighter sources having a higher multiplicity rate. 
This result is also consistent with multiplicity studies 
of 5'870-selected DSFGs by Ivison et al. (2007), Smolcic 
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Figure 5. Cumulative distribution functions showing the effect of 
lensing on the inferred flux densities of the Herschel-ALMA sample. 
The median flux density in the Herschel-AhMA sample drops by 
a factor of 1.3 when lensing is taken into account. For comparison 
we also show the flux density distribution from ALESS (pink line), 
with the completeness limit of the LESS survey indicated by a 
dashed pink line. 

et al. (2012), and Barger et al. (2012), who use VLA, 
PdBI/1.3mm, and SMA/870/im imaging to determine 
rates of 18%, 22%, and 40%, respectively. 

One useful way to characterize multiplicity is with 
a comparison of the total 870/im flux density, 5'totaL 
with the individual component 870 fim flux density, 
5'component- Figure 6 shows these values for our Herschel- 
ALMA sample and compares to ALESS. Lensing has a 
significant impact on the apparent flux densities of many 
objects in our ALMA sample, so we are careful to show 
only intrinsic flux densities in this diagram. This dia¬ 
gram reflects the known result that the multipicity rate 
in ALESS rises and the average fractional contribution 
per component decreases with increasing 5'totai (Hodge 
et al. 2013). A simple extrapolation of this phenomenon 
to the flux density regime probed by our Herschel-ALMA 
sample would have suggested a very high multiplicity rate 
and a very low average fractional contribution per com¬ 
ponent. The multiplicity rate in our sample is indeed 
higher, but we find that the average fractional contribu¬ 
tion per component hovers around 0.4 for essentially the 
full range in our sample. This is a reflection of the fact 
that the brightest Herschel DSEGs comprise 1-3 ALMA 
components, not 5-10 ALMA components as might have 
been expected from a naive extrapolation of the ALESS 
results. 

4.3. Spatial Distribution of Multiple Sources 

We can dig further into our ALMA data by exploring 
the average number of ALMA sources per annular area 
(dN/dA) as a function of how far they are from each 
other. Eigure 7 shows the results of this analysis for 
both our Herschel-ALMA sample and ALESS. We formu¬ 
late the separation as an angular distance between each 
ALMA source (using the lensing-corrected data) and the 
centroid of all of the ALMA sources for that Herschel 



Figure 6. Comparison of the total 870/rm flux density, ^totab 
with the individual component 870 p,m flux density, S'component 
(both of these are after accounting for lensing). Objects falling 
along the gray dashed line are single component systems (i.e., 
'8totai = 8'component)- The solid lines trace the average ratio of 
component to total flux for a given total flux. Our sample of Her¬ 
schel DSFGs {Herschel-AIM A sample, green squares) has a higher 
multiplicity and a lower average factional contribution per compo¬ 
nent than the ALESS sample (pink diamonds), but not as low as 
would be expected from a simple extrapolation of the trend in the 
ALESS data alone. 

DSEG. This is different from the pairwise separation dis¬ 
tance estimator used by Hodge et al. (2013) that be¬ 
comes ill-defined when there are more than two ALMA 
counterparts (as is often the case in our Herschel-ALMA 
sample). Eigure 7 shows dNjdA values for ALESS that 
have been re-computed using our method. We also show 
the median and la range found from simulated datasets 
for both ALESS and our Herschel-ALMA sample. The 
simulated datasets consist of 200 runs of DSFGs with 
the same flux density and multiplicity as the observed 
datasets (both the ALESS sample and our ALMA sam¬ 
ple), but placed randomly within the primary beam 
FWHM. We also show predictions from simulations by 
HB13 (see below for details). 

We recover the result from Hodge et al. (2013) that 
the ALESS DSFGs are consistent with a uniformly dis¬ 
tributed population. Interestingly, however, there is a 
dramatic rise in dN/dA for angular separations less than 
2" in our Herschel-ALMA sample. Indeed, for an an¬ 
gular separation of 075, we find an excess in dN/dA 
by a factor of 10 compared to a random, uniformly 
distributed population. This excess persists (although 
at significantly lower amplitude) even when the quality 
of our ALMA observations are degraded to match the 
typical sensitivity, spatial resolution, and uv coverage of 
ALESS (as represented by observations of ALESS 122). 
The persistence of the excess suggests that it is an intrin¬ 
sic property of the sample; i.e., that only the brightest 
DSFGs show an excess of sources on small separation 
scales (with the caveat that we cannot rule out the pos¬ 
sibility of at least part of the excess being due to strong 
lensing from optically dark lenses). 

An excess of sources with small separations from each 
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Figure 7. Left: Spatial distribution of sources with multiple counterparts found in our Herschel-ALMA sample (green squares), in ALESS 
(pink diamonds) and mock catalogs from HB13 and C15 (blue plus signs and orange crosses, respectively). Sources identified in our 
Herschel-ALMA sample lie much closer to each other than they do in either ALESS or the HB13 simulations. Right: Number of ALMA 
sources per annular area as a function of angular separation from the ALMA centroid. Symbols and colors are as in left panel. We also show 
how our Herschel-ALMA sample would appear if it had been observed with ALESS resolution and sensitivity (light green squares). The 
range of separations that would be seen if sources were randomly distributed within the ALMA field of view is also shown (dark dashed line 
and hatched region). The Herschel-ALMA DSFGs show a significantly stronger excess on angular separation scales < 2" compared to both 
ALESS and the HB13 simulations, even when taking into account the difference in sensitivity and spatial resolution between our ALMA 
observations and those of ALESS. The simulations from C15 show better agreement with the data, likely due to the more sophisticated 
treatment of blending compared to HB13. 


other could be an indication of interacting or merging 
systems. However, it is also possible that the sources 
are merely unrelated galaxies that appear blended due 
to projection effects (Hayward et al. 2013a; Cowley et al. 
2015; Munoz Arancibia et al. 2015). Spatially resolved 
spectroscopy is necessary to answer this question defini¬ 
tively, but is not currently available. Instead, to inves¬ 
tigate these possibilities further, we make use of mock 
catalogs of DSFGs that are based on numerical simula¬ 
tions and presented by HB13 and C15. 

We begin with the HB13 simulations, summarizing the 
methodology used to generate the mock catalogs here 
and refering the reader to HB13 for full details. 

Halo catalogs are generated from the Bolshoi dark 
matter-only cosmological simulation (Klypin et al. 2011) 
using the ROCKSTAR halo finder (Behroozi et al. 2013b,c). 
Catalogs of subhalos are created from eight randomly 
chosen lightcones, each with an area of 1.4° x 1.4°. 
Galaxy properties such as stellar mass and SFR are 
assigned to the subhalos using the abundance match¬ 
ing method of Behroozi et al. (2013a). Dust masses 
are assigned using an empirically determined redshift- 
dependent mass-metallicity relation and an assumed 
dust-to-metal density ratio of 0.4 (see Hayward et al. 
2013b for details). Finally, submm flux densities are in¬ 
terpolated from the SFRs and dust masses using a fitting 
function that is based on the results of dust radiative 
transfer calculations performed on hydro dynamical sim¬ 
ulations of isolated and interacting galaxies (Hayward 
et al. 2011, 2012b, 2013b). 

A blended source is deflned as any galaxy in the mock 
catalogs above a threshold flux density (5'thresh) that 
has at least one neighbor within a projected angular 
distance ^neighbor- To obtain a direct comparison with 


our Herschel-ALMK sample, we use 5'thresh = 1.0 mJy 
(corresponding to the 5cr limit of the ALMA data) and 
^neighbor = 40" (reflecting the size of the Herschel beam 
at 500 /im). We use the known positions in the mock 
catalogs for all blended sources and compute centroid 
and separations for every blended source using the same 
methodology as we applied to our Herschel-KhMA sam¬ 
ple and to ALESS. 

The dN/dA values found in the mock catalogs are 
shown by the thick blue line and plus signs in Figure 7. 
There is a significant increase in dNjdA on separations 
smaller than 075, but the amplitude of the increase 
is much lower than is apparent in our Herschel-ALMA 
sample. 

The HB13 model does not include SFR enhancements 
induced by starbursts (see Section 4.5 of Hayward et al. 
2013a for a detailed discussion of this limitation). To ex¬ 
plore whether interaction-induced starbursts are the ori¬ 
gin of the excess at small angular separations observed in 
our Herschel-ALMA sample, we analyzed modified ver¬ 
sions of the HB13 model that include a crude treatment 
of interaction-induced SFR enhancements (Miller et al. 
2015). Mock galaxies with one or more neighbors within 
a physical distance of 5kpc and with a stellar mass be¬ 
tween one-third and three times that of the galaxy under 
consideration (i.e. a ‘major merger’) had their SFRs in¬ 
creased by a factor of two. For distances smaller than 
1 kpc, the imposed increase was a factor of ten. Because 
these SFR enhancements are greater than suggested by 
simulations (e.g. Cox et al. 2008; Hayward et al. 2011, 
2014; Torrey et al. 2012) or observations of local galaxy 
pairs (e.g. Scudder et al. 2012; Patton et al. 2013), we 
consider this test to provide an upper limit on the possi¬ 
ble effect of interactions on blended sources in the HB13 
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model, although the incompleteness of the Behroozi et al. 
(2013c) catalogs for mergers with small separations could 
cause some interacting systems to be missed. We find an 
insignificant effect on the values of dN/dA when using 
the merger-induced model as described above. The main 
reason for this is that only two sources had their SFRs 
boosted by a factor of ten, and ^150 experienced a fac¬ 
tor of two increase. In the HB13 model, a factor of two 
increase in SFR corresponds to only a 30% increase in 
Ss 70 i so it is perhaps unsurprising that the weak boosts 
in SFR cause little change in dN/dA. 

Experiments with stronger interaction-induced SFR 
enhancements showed that very high enhancements (e.g. 
a factor of 10 for separation of 5-15 kpc and 100 for sepa¬ 
ration of < 5 kpc) in major mergers are required to match 
the observed excess in dN/dA on small separations. In¬ 
corporating starbursts induced by minor-merger could 
possibly reduce the required SFR enhancements. The 
tension between the model prediction and observations 
may also indicate that a more sophisticated treatment of 
blending is necessary. 

To explore this possibility, we investigate mock cata¬ 
logs based on the methodology presented in C15. Here, 
we give a brief summary and refer the interested reader 
to C15 for full details. A new version of the galform 
(e.g. Cole et al. 2000, Lacey et al. in preparation) semi- 
analytic model of hierarchical galaxy formation is used 
to populate halo merger trees (e.g. Parkinson et al. 2008; 
Jiang et al. 2014) derived from a Millennium style N- 
body dark matter only simulation (Springel et al. 2005; 
Guo et al. 2013) with WMAP7 cosmology (Komatsu 
et al. 2011). A sub-mm flux for each galaxy is calculated 
using a self-consistent model based on radiative trans¬ 
fer and energy balance arguments. Dust is assumed to 
exist in two components, dense molecular clouds and a 
diffuse ISM. Energy absorbed from stellar radiation by 
each dust component is calculated by solving the equa¬ 
tions of radiative transfer in an assumed dust geometry. 
The dust is then assumed to emit radiation as a modified 
blackbody. 

Three randomly orientated 20 deg^ lightcone cata¬ 
logues are generated using the method described in Mer- 
son et al. (2013). We choose as the lower flux limit for 
inclusion of simulated galaxies into our lightcone cata¬ 
logue S'sooAim >0.1 mJy, as this is the limit at which we 
recover 90% of the extragalactic background light (EBL) 
as predicted by the model (122 Jy deg“^). This is in 
excellent agreement with observations from the COBE 
satellite (e.g., Puget et al. 1996), and thus ensures a re¬ 
alistic 500 /rm background in the mock images. 

Mock imaging is created by binning the lightcone 
galaxies onto a pixelated grid which is then convolved 
with a 36" EWHM Gaussian (corresponding to the iJer- 
schel/SFIRE beam at 500//m). The image is then con¬ 
strained to have a zero mean by the subtraction of a 
uniform background. No instrumental noise is added, 
nor are any further hltering procedures applied to the 
mock image. Eor the purposes of source identihcation, 
this procedure is repeated at 250 jam. Eor this, we adjust 
the EWHM of the Gaussian PSE to 18" and change the 
lower flux limit of inclusion into our lightcone to ensure 
90% of the predicted EBL is recovered at this wavelength. 

Source positions are selected as maxima in the mock 


250 /rm image, with the position of the source being 
recorded as the center the maximal pixel for simplicity. 
To mimic ‘deblended’ Herschel photometry we record the 
value of the pixel located at the position of the 250 /im 
maxima in the 500 /im images. We select ah Herschel 
sources satisfying S'soo/xm > 50mJy and z > 1 to identify 
galaxies from our lightcone catalogs within a 9" radius of 
the source position, modelling the ALMA primary beam 
prohle as a Gaussian with an 18" EWHM and a maxi¬ 
mum sensitivity of 1 mJy. 

The dN/dA values derived from the G15 mock catalogs 
are shown by the thick orange line and crosses in Eig- 
ure 7. Here, the amplitude of the increase in dN/dA on 
separations smaller than ^ 2" mimics the trend seen in 
the data. However, there is a dehcit of multiple systems 
with separations of 075 or less compared to the Herschel- 
ALMA sample. This result suggests that a sophisticated 
treatment of blending yields better agreement between 
simulations and observations but the simulations still 
under-predict the number of multiple systems with small 
separations. 

Euture work on the theoretical side should seek to de¬ 
termine if the application of the G15 blending algorithm 
to the HB13 simulations yields similarly better agree¬ 
ment with the data. On the observational side, it is crit¬ 
ical to establish whether Herschel sources with multiple 
ALMA counterparts are physically related by measuring 
spectroscopic redshifts to individual counterparts. Eor- 
tunately, this is a viable project today with the VLA and 
ALMA. 

5. IMPLICATIONS FOR THE BRIGHT END OF THE DSFG 
LUMINOSITY FUNCTION 

The distribution of magnihcation factors for sources 
found in wide-held surveys with the brightest apparent 
hux densities are highly sensitive to the shape of the 
intrinsic number counts at the bright end. In this sec¬ 
tion, we combine our ALMA and SMA measurements of 
magnihcation factors to investigate this as it pertains to 
DSEGs. 

5.1. Statistical Predictions for jasjo 

Our methodology follows the procedures outlined in 
previous ehorts to predict magnihcation factors for DS¬ 
EGs with a given apparent hux density (chiehy, Lima 
et al. 2010; Wardlow et al. 2013; Eialkov & Loeb 2015). 
We summarize the essential elements here and highlight 
signihcant diherences where appropriate. 

The key components of the model are the mass den¬ 
sity prohle of the lenses, piens the number density of 
lenses as a function of mass and redshift, niens z), the 
redshift distribution of the sources, dNsource/dz, and the 
intrinsic number counts of the sources, dUsource/dSgYo- 
The latter component is the least certain and also has 
the strongest impact on the predicted apparent luminos¬ 
ity function. Eor these reasons, we hx ah components 
of the model except the shape of the intrinsic number 
counts. Our goal is to take luminosity functions that 
can successfully ht observed faint DSEG number counts 
(Karim et al. 2013) and test whether they lead to pre¬ 
dicted magnihcation factors consistent with our ALMA 
and SMA observations. 

To describe piens (^)5 we use a superposition of a singu¬ 
lar isothermal sphere (SIS) and a Navarro-Erenk-White 
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(NFW) profile (Navarro et al. 1997) that is truncated at 
the virial radius. The NFW profile describes the out¬ 
skirts of dark matter halos better (Mandelbaum et al. 
2005), while the SIS profile is preferred on smaller scales 
because it correctly fits the observed fiat rotational 
curves in galaxies (Kochanek 1994). We make sure that 
the resulting probability density of lensing, P (/r), is nor¬ 
malized to unity. To describe niens z), we generate 
the abundance of halos at each mass and redshift us¬ 
ing the Sheth & Tormen (1999) formalism. To describe 
dNsource/d^^ wc adopt the following redshift distribution 
which is based on photometric redshifts of optical coun¬ 
terparts to ALMA sources identified in ALESS (Simpson 
et al. 2014): 


dN/dz oc 


azCTzV^ 


exp 


-[In(a^) -ln(l + 2 ;^)]^ 
2(jlaz 


( 2 ) 

where = 1 -\- z^ = 1.5 (to reflect the relatively blue 
SPIRE colors of the sample), and (Jz = 0.2. Alterna¬ 
tive values for z^ and yield second-order perturba¬ 
tions which are not significant at the level of our current 
analysis. 

We explore two intrinsic number counts that are in¬ 
tended to bracket the plausible range of values for DS- 
EGs based on two interferometric surveys. One is the 
number counts measured in ALESS (Karim et al. 2013), 
and the other is from interferometric follow-up of the 
first AzTEC survey in COSMOS (Scott et al. 2008) using 
the SMA (Younger et al. 2007, 2009) and PdBI (Mietti- 
nen et al. 2015). These interferometric observations have 
shown that all of the sources in their surveys are either 
unlensed or lensed by magnification factors < 2 (a simi¬ 
lar result is found based on ALMA imaging of 52 DSEGs 
in the Ultra Deep Survey; Simpson et al. 2015a). This 
is why the ALESS and COSMOS luminosity functions 
represent a plausible range of intrinsic number counts for 
DSEGs. These number counts are shown in the left panel 
of Eigure 8. Interferometric follow-up data in COSMOS 
(Smolcic et al. 2012) and GOODS-N (Barger et al. 2012) 
are published, but unknown completeness corrections in 
the single-dish surveys on which these follow-up datasets 
are based precludes their use here. 

In detail, we use a broken power-law of the form 


distribution function. 


Table 4 

Parameters of DSFG luminosity functions tested in this paper. 


Luminosity Function 

Af* 

(deg-2) 

(mJy) 

/9i 

P2 

ALESS broken power-law 

20 

8 

2 

6.9 

GOSMOS broken power-law 

20 

15 

2 

6.9 


The fundamental measurement provided by the spa¬ 
tially resolved ALMA and SMA imaging and associated 
lens models is the magnification factor of a source with 
a given apparent S'gTO- We use the combined sample to 
compute the average magnification as a function of S'gTO 
from the data: (// 87 o)- The same quantity can also be 
directly computed from our model as 

pOO 

(M87o) = / fiP{fi\Ss 7 o)dfi, (4) 

Jo 

where the probability for lensing with magnification ji 
given the apparent fiux is: 


and 


P{fl\Ss7o) 


1 p(/i) dn / Ss 7 o \ 
N dS',,, W / ’ 


j /i dSgjQ \ J 


(5) 

(6) 


Here dn/dSs 7 o is the observed number counts and 
dn/dSgjQ is the intrinsic number counts. 

As part of the lens models, the ALMA and SMA imag¬ 
ing also provide the probability that a source with a given 
apparent 5'870 experiences a magnification above some 
threshold value, /imin- P {/Ji > Mmin)- It is therefore of 
interest to make a similar prediction from our model. We 
use the following to do this: 


P{/Jj ^ /^min) 


/y p{ii\Sobs) 
j^piiiiSobsY 


(7) 


Table 4 provides values for the parameters of the bro¬ 
ken power-law for the ALESS and GOSMOS number 
counts. The data and corresponding number counts are 
shown in the left panel of Eigure 8. 

The product of the model is the lensing optical depth 
for a given lensing galaxy and source galaxy, /^. The 
lensing probability with magnification larger than /i is 
then calculated via P{> fi) = 1 — exp(—/^) and the 
differential probability distribution is P(/i) = —dP{> 
[i)jd[i. The sum over the distribution of source redshifts 
and lens masses and redshifts yields the total probability 


5.2. Comparing Models with Data 

The middle panel of Eigure 8 shows a direct compari¬ 
son of the measured /i 87 o values as a function of apparent 
Ss7o for the Herschel-AljM.A and Herschel-^yiA samples. 
We also show a running average of the combined sample 
(considering only yU 870 > 2.0 objects) to serve as a di¬ 
rect comparison to our theoretical models. We compute 
this by interpolating the observed /r 87 o and 5'870 onto a 
fine grid using the SciPY griddata package and then 
smoothing the resulting grid using a Gaussian filter in 
the Scipy gaussian_filter package. Also shown in this 
diagram are model predictions for the average magnifi¬ 
cation as a function of S^ 7 q^ (^ 870)5 assuming the two 
intrinsic number counts for DSEGs described in Table 4. 

Both models predict higher (// 87 o) than are seen in 
the data by factors of 5-10. However, the dispersion 
in the predicted (/i 87 o) values for both intrinsic number 
counts rises smoothly from ^ 2 at 5'870 = 15mJy to 
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Figure 8. Left: Observed number counts from interferometer follow-up of mm sources in COSMOS (black circles; Younger et al. 2007, 
2009; Miettinen et al. 2015) and from ALESS (pink diamonds; Karim et al. 2013). The magneta and blue lines represent the range of 
plausible intrinsic number counts for DSFGs given these two datasets. Middle: Magnification factors at 870/rm as a function of apparent 
Ss 7 o for every source in our ALMA (green squares) and SMA (red circles) samples. The black line represents a running average of the 
magnification measurements when sources with p-sro > 2.0 are considered (grey shaded region highlights the Icr uncertainty), and is 
truncated at the minimum fiux density at which we have observed figro > 2.0. Colored lines show the average magnification factor, (/isro), 
predicted by our models for the two intrinsic number counts shown on the left. The model number counts predict higher {/J^sro) than are 
seen in the data. Right: Probability that a source with a given Sgro experiences /isTO > 2.0. The black line (grey shaded region) shows a 
logistic regression fit to the SMA and ALMA data (grey shaded region highlights the Icr uncertainty). Colored lines are the same as in the 
middle and left panels. The intrinsic number counts that provide a good fit to the COSMOS data predict too many unlensed or weakly 
lensed sources with intrinsic fiux densities of Sg^Q ~ 50 mjy. 


« 8 at 5'87 o = 100 mJy, so this difference is not statis¬ 
tically significant. Furthermore, there is reason to believe 
that the data may be biased against high magnification 
factor measurements. In both the Herschel-ALMA and 
Herschel-SMA samples, the spatial resolution is ~ 0'.'5. 
This is nearly always sufficient to resolve the images of 
the lensed galaxy, but it is not always sufficient to re¬ 
solve the images themselves. Therefore, it may be the 
case that the lens models over-predict the intrinsic sizes 
of the lensed galaxies and hence under-predict the mag¬ 
nification factors. For example, the average half-light 
sizes of unlensed DSFGs have been found recently to be 
^ 0.8 kpc (Ikarashi et al. 2014) and « 1.2 kpc (Simpson 
et al. 2015a). In contrast, we reported a median half- 
light radius of 1.53kpc in Bussmann et al. (2013). Lens 
models from higher resolution data with ALMA suggest 
that magnification factors could increase by a factor of 
^ 1.5 — 2 (Rybak et al. 2015; Tamura et al. 2015; Dye 
et al. 2015). Therefore, it is plausible that both of the 
intrinsic luminosity functions tested here provide statis¬ 
tically consistent fits to the data. 

A related but distinct test of the intrinsic number 
counts for DSFGs comes from the probability of a given 
source experiencing a magnification above some thresh¬ 
old value, /imin- Unlike the case with the average magni¬ 
fication factor measurements, our ALMA and SMA data 
should provide a reliable estimate of this quantity. The 
results of this are shown in the right panel of Figure 8. 
For clarity of presentation, we show one choice of /imin- 
Rmin = 2.0. The shape of the curves varies with /imin, 
but the overall results are qualitatively the same. The 
models we consider are the same as those used in the 
left panel of Figure 8. Instead of computing a running 
average of the data, we show a logistic regression fit to 
the data (obtained with the Scikit-learn package; Pe- 
dregosa et al. 2011). 

Both models tested in this paper exhibit a sharp tran¬ 
sition from low probability to high probability of being 
lensed, consistent with the data. However, there are sig¬ 


nificant differences in where this transition fiux density, 
S'trans, occurs — i.c., whcrc P(/i > R 87 o) >2. In the 
data, the logistic regression fit yields Ptrans = 24 ± 3 mJy 
(error accounts only for statistical uncertainty in SgTo 
and /i 87 o), whereas the models based on the ALESS 
and COSMOS number counts yield Ptrans = 37 mJy 
5'trans = 69 mJy, respectively. 

This analysis highlights the difficulty encountered with 
the luminosity function based on the COSMOS data: un¬ 
lensed sources with S '870 > 50 mJy are over-predicted and 
lensed sources with SgTo < 50mJy are under-predicted. 
If the ALESS number counts continue to be supported 
by the evidence as additional data are obtained (e.g.; 
Simpson et al. 2015b), the implications are significant. 
We should then expect to find ~ 3 sources satisfying 
SgTo > 10 mJy in a 1 deg^ survey. This is about a fac¬ 
tor of 7 lower than typical measurements from single¬ 
dish, broad-beam studies (e.g., WeiB et al. 2009). This 
suggests that very luminous galaxies such as GN20 and 
HELS3 may be more rare than previously thought. 

6. CONCLUSIONS 

We present ALMA 870 /im 0745 imaging of 29 Herschel 
DSEGs selected from 55deg^ of HerMES. The Herschel 
sources have S^qq = 52 — 134 mJy, placing them in a 
unique phase space between the brightest sources found 
by Herschel and those found in ground-based surveys at 
sub-mm wavelengths that include more typical, fainter 
galaxies. Our ALMA observations reveal 62 sources 
down to the 5cr limit (a ~ 0.3mJy, typically). We make 
use of optical and near-IR imaging to assess the distri¬ 
bution of intervening galaxies along the line of sight. We 
introduce a new, publicly available software called UVM- 
CMCFIT and use it to develop lens models for all ALMA 
sources with nearby foreground galaxy. Our results from 
this effort are summarized as follows: 

1. 36/62 ALMA sources experience significant ampli¬ 
fication from a nearby foreground galaxy that is 
comparable to or greater than the absolute call- 
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bration uncertainty (i.e., /igro > l-l)- The median 
amplification in the subset that experiences lens- 
ing is /i 87 o = 1-b. Only 6 sources show morphology 
typical of strong gravitational leasing and could be 
identified as lenses from the ALMA imaging alone. 
A multi-wavelength approach is critical to identify¬ 
ing structure along the line of sight and determin¬ 
ing an unbiased measurement of the flux densities 
in our sample. 

2. 20/29 Herschel DSFGs break down into multiple 
ALMA counterparts. Of the 9 isolated systems, 5 
are strongly leased by factors of 5-30 (HECDFS12 
is a non-isolated system with a strongly leased 
source in it). After correcting for amplification, 
the brightest source in the sample has Ssto = 
17.5 ± 0.4 mJy. There is a weak trend towards even 
higher multiplicity at the highest total Ssto flux 
densities. 

3. When a Herschel source comprises multiple ALMA 
counterparts, these counterparts are typically lo¬ 
cated within 2" of each other. Their separations 
are significantly smaller than ALMA counterparts 
to ALESS sources as well as simulated sources from 
HB13 and C15, although the improved treatment 
of blending by the latter yields superior agreement 
with the data. This conclusion remains true even 
when we degrade our ALMA observations to match 
the spatial resolution, sensitivity, and uv coverage 
of the ALESS observations. 

4. Intrinsic number counts for DSEGs with a form 
that matches observations in GOSMOS (Miettinen 
et al. 2015) under-predict the number of lensed 
sources with apparent Ssto > lOmJy. Number 
counts based on ALESS observations provide a bet¬ 
ter match to our magnification measurements. The 
interpretation of these results is complicated by the 
fact that our sample is likely biased towards blends 
of multiple sources within a Herschel beam. Our 
primary goal is to draw attention to this analysis 
method as a means to test number counts of DS¬ 
EGs using wide-held Herschel data. 

If the ALESS number counts continue to provide the 
best predicted magnihcation factors in larger samples 
with weaker biases, this suggests that galaxies with in¬ 
trinsic hux densities above ^870 ^10 mJy are extremely 
rare. One possible explanation for their rarity is that 
they are simply the tip of the mass function among star- 
bursts. An alternative is that they represent a very short 
phase in galaxy evolution. It is interesting to note that 
consistent with this idea is the high multiplicity rate in 
our ALMA sample as well as the small projected sepa¬ 
rations between multiple ALMA counterparts. The in¬ 
ability of numerical simulations to reproduce the small 
projected separations seen in the data might highlight 
a productive path forward to improve our theoretical 
understanding of the enhancement in star-formation by 
interactions and mergers of galaxies which are already 
forming stars at a very high rate. 

In the future, higher spatial resolution imaging is 
needed to investigate the morphologies of individual 
ALMA sources. Tidal tails, multiple nuclei, and other 


signs of mergers and interactions should become evident 
at 071 resolution. In addition, molecular spectroscopy 
will be critical to determine distances to and dynam¬ 
ics of individual ALMA sources and hence characterize 
what fraction of Herschel sources are actually physically 
associated with each other (and not simply a result of 
projection effects along the line of sight). Einally, deeper 
optical or near-IR imaging is needed to search for dark 
lenses that may have been missed in our Gemini-S or 
HST imaging. 

We thank the referee for useful comments that helped 
improve the clarity of the manuscript. This paper makes 
use of the following ALMA data: ADS/JAO.ALMAt^ 
2011.0.00539.S. ALMA is a partnership of ESO (rep¬ 
resenting its member states), NSE (USA) and NINS 
(Japan), together with NRG (Ganada) and NSG and 
ASIA A (Taiwan), in cooperation with the Republic of 
Ghile. The Joint ALMA Observatory is operated by 
ESO, AUI/NRAO and NAOJ. 

The results described in this paper are based on obser¬ 
vations obtained with Herschel^ an ESA space observa¬ 
tory with science instruments provided by European-led 
Principal Investigator consortia and with important par¬ 
ticipation from NASA. 

This research has made use of data from the Her- 
MES project (http://hermes.sussex.ac.uk/). HerMES is 
a Herschel Key Programme utilizing Guaranteed Time 
from the SPIRE instrument team, ESAG scientists, and 
a mission scientist. HerMES is described in Oliver et al. 
(2012). The HerMES data presented in this paper will 
be released through the Herschel Database in Marseille 
(HeDaM^i). 

SPIRE has been developed by a consortium of in¬ 
stitutes led by Gardiff Univ. (UK) and including: 
Univ. Lethbridge (Ganada); NAOG (Ghina); GEA, LAM 
(Erance); lESI, Univ. Padua (Italy); lAG (Spain); Stock¬ 
holm Observatory (Sweden); Imperial Gollege London, 
RAL, UGL-MSSL, UKATG, Univ. Sussex (UK); and 
Galtech, JPL, NHSG, Univ. Golorado (USA). This devel¬ 
opment has been supported by national funding agencies: 
GSA (Ganada); NAOG (Ghina); GEA, GNES, GNRS 
(Erance); ASI (Italy); MGINN (Spain); SNSB (Sweden); 
STEG, UKSA (UK); and NASA (USA). 

The VIDEO imaging data are based on data products 
from observations made with ESO Telescopes at the La 
Silla or Paranal Observatories under ESO programme ID 
179.A-2006. 

This paper is partly based on observations obtained at 
the Gemini Observatory, which is operated by the Asso¬ 
ciation of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSE on behalf 
of the Gemini partnership: the National Science Eoun- 
dation (United States), the National Research Gouncil 
(Ganada), GONIGYT (Ghile), the Australian Research 
Council (Australia), Ministerio da Ciencia, Tecnologia e 
Inovacao (Brazil) and Ministerio de Ciencia, Tecnologfa 
e Innovacion Productiva (Argentina). 
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APPENDIX 

HERSCHEL/SPIRE PHOTOMETRY 

We present in Table 5 SPIRE photometry {S 250 , 'S'sso, and S'soo) for each Herschel source in our sample using as 
priors the ALMA and Spitzer 24 /im counterpart positions. Many of the Herschel sources in our sample have multiple 
24/rm counterparts close enough to the ALMA counterparts that they can make a significant contribution to the 
SPIRE flux at the position of the ALMA counterparts. For this reason, it is critical to include the 24 /rm data when 
estimating Herschel/SFIRE photometry. In cases where the Spitzer 24/rm and ALMA counterparts spatially overlap 
(defined here as having a separation smaller than 2'.'5), we exclude the Spitzer counterpart from the calculations. 

Note that in the case of HELAISS02, we exclude a 24/rm source associated with the leasing galaxy, despite being 
more than 2'.'5 away from the nearest ALMA counterpart. This is because the ALMA counterparts are arranged so 
that they surround the 24/im source completely. The extended emission seen in the SPIRE maps is therefore most 
likely attributable to the ALMA counterparts rather than the 24 /im source. 

One of the key results from this analysis is that the use of StarFinder to deconvolve the SPIRE beam will lead to 
the removal of a portion of the targets with extended emission from a flux-limited sample. This is by design, as the 
goal outlined in Wardlow et al. (2013) was to develop the purest sample of lens candidates as possible from HerMES 
data. Blends were considered by these authors to be unlikely to be lensed. This method is effective for lenses where 
the Einstein radius is smaller than ^ 2". However, it selects against deeper potential wells, such as those of groups 
or clusters (e.g., HLockOl Gavazzi et al. 2011), that produce images separated by scales that comprise a significant 
fraction of the SPIRE PSF. 

It is worth emphasizing that blends likely constitute an interesting path of study for future work, as they potentially 
represent proto-groups or proto-clusters during a particular active stage of their evolution. They are also prime 
examples of systems of sources that are poorly reproduced in simulations, as evidenced by the investigations in 
Section 4.3 in this paper. Table 5 shows that one possible means of selecting candidate blends is by comparing 
the SUSSEXtractor and StarFinder flux densities, particularly at 500/im. A large difference between these two 
measurements likely indicates multiple components separated by scales that are a significant fraction of the SPIRE 
beam (e.g., HELAISS02, HXMM20, and HGOSMOS02). 

Table 5 

Gompilation of Herschel/SPIRE flux density measurements. For each 
Herschel source, we give the flducial flux densities (denoted in table as 
“Fiducial”), the initial flux densities obtained using SUSSEXtractor that 
were then used to generate the target list for the ALMA observations 
(“SUSSEXtractor”), and the flux densities obtained subsequently using the 
more sophisticated deblending algorithm from StarFinder that were then 
used to construct the list of lens candidates presented in Wardlow et al. 

(2013) (“StarFinder”). In all cases, uncertainties are comparable to the 
uncertainties given in Table 1. 


Short name 

‘S '250 

(mjy) 

Fiducial 

^'350 

(mJy) 

‘S'soo 

(mJy) 

S 250 

(mJy) 

StarFinder 

‘S '350 

(mJy) 

‘S'soo 

(mJy) 

SUSSEXtractor 

S 250 ‘S '350 S'soo 

(mJy) (mJy) (mJy) 

HELAISS02 

115 

124 

108 

114 

101 

76 

105 

128 

103 

HXMM02 

78 

122 

99 

92 

122 

113 

101 

147 

141 

HXMM31 

102 

94 

65 

128 

112 

73 

129 

116 

80 

HXMM29 

97 

102 

78 

89 

83 

56 

100 

107 

80 

HXMM07 

89 

107 

85 

91 

104 

86 

92 

104 

83 

HXMM20 

72 

85 

66 

85 

79 

67 

80 

96 

88 
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Fiducial StarFinder SUSSEXtractor 


Short name 

‘S'250 

(mJy) 

‘S'350 

(mJy) 

‘S'500 

(mJy) 

‘S'250 

(mJy) 

‘S'350 

(mJy) 

‘S'500 

(mJy) 

‘S'250 

(mJy) 

‘S'350 

(mJy) 

‘S'500 

(mJy) 

HXMMOl 

179 

188 

134 

180 

192 

132 

178 

195 

137 

HXMM04 

162 

157 

125 

144 

137 

94 

173 

174 

127 

HXMM09 

129 

118 

85 

120 

115 

84 

135 

113 

95 

HXMM03 

114 

134 

116 

121 

132 

110 

118 

137 

118 

HXMMll 

101 

104 

73 

107 

108 

81 

105 

121 

94 

HXMM23 

128 

105 

68 

137 

108 

57 

132 

118 

88 

HXMM22 

101 

85 

61 

97 

82 

62 

147 

128 

89 

HXMM05 

103 

118 

97 

106 

119 

92 

103 

115 

101 

HXMM30 

86 

97 

75 

90 

100 

75 

93 

105 

80 

HXMM12 

98 

106 

82 

102 

110 

81 

107 

115 

89 

HECDFS12 

61 

82 

81 

28 

84 

85 

68 

92 

100 

HECDFS04 

56 

61 

55 

73 

86 

85 

65 

87 

96 

HECDFS13 

95 

89 

63 

96 

90 

63 

88 

85 

51 

HADFS04 

74 

93 

84 

76 

90 

72 

71 

95 

87 

HADFS02 

102 

97 

81 

110 

102 

87 

103 

100 

79 

HADFSll 

19 

39 

52 

57 

78 

75 

57 

87 

97 

HADFSIO 

47 

58 

58 

96 

86 

57 

121 

114 

76 

HADFSOl 

76 

100 

94 

80 

103 

93 

72 

108 

87 

HADFS09 

98 

102 

72 

115 

61 

24 

112 

117 

86 

HADFS08 

142 

133 

90 

88 

81 

50 

126 

130 

102 

HADFS03 

119 

102 

63 

138 

114 

73 

134 

124 

86 

HCOSMOS02 

70 

85 

71 

71 

64 

41 

82 

99 

89 

HCOSMOSOl 

86 

96 

71 

91 

100 

74 

89 

99 

73 



